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Foreword 

This Handbook was prepared by Northrop Space Laboratories, 
Huntsville, Alabama, for the Dynamics and Analysis Section of 
the Aero-Astrodynamics Laboratory of the George C. Marshall 
Space Flight Center. Its contents represent a partial fulfill- 
ment of the project since additional topics are forthcoming 
which, upon completion, will be integrated into the Handbook. 
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Abstract 


3 ? 3 3 ^ ^ 


Design information related to the effects of propellant sloshing 
is presented for use in both cont^ 1 °nd structural problems. Both 
analytical and experimental results are given and all pertinent material 
is referenced. Graphs have been included, whenever possible, to expedite 
preliminary design calculations. The areas covered are: (1) linearized 

fluid theory, (2) equivalent mechanical model theory, (3) results of 
analytical studies of liquid oscillations in variously shaped containers 
when subjected to different types of excitation, i.e., boundary conditions, 
fluid velocity potentials, natural frequencies, liquid force and moment 
resultants and equivalent mechanical models, and (4) results of both 
analytical and experimental studies concerned with propellant slosh 
suppression, with particular emphasis on fixed-ring baffles. 
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NOTATION 
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F 

F(s) 

F(t) 
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K 

H 

n 

h 

h 
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Radius of circular cylindrical tank 

Outer radius of annular or annular-sector tank 

Length of rectangular tank in x-direction 

Inner radius of ring baffle 

Inner radius of annular or annular-sector tank 
Local drag coefficient 

Radial clearance between tank wall and ring baffle 

Effective damping coefficient 

Damping coefficient of disc 

Damping coefficient of n^ slosh mass 

Center of gravity of undisturbed liquid 

Baffle spacing 
Tank diameter 

Tank diameter 

Depth of baffle below undisturbed free surface 

Depth of baffle below undisturbed free surface 

Distance of baffle above the undisturbed free surface, 
D-d 

Force in q. -direction 
1 

3 

Dimensionless liquid force, F/pga 

Laplace transform of the force F 

Arbitrary function of time 

Acceleration due to gravity 

th 

Damping factor of n liquid mode 

z-coordinate for point of attachment of n^ pendulum 

Depth of liquid measured from tank bottom 
th 

z-coordinate of n slosh mass 
Mass moment of inertia of fixed mass 
vi 



NOTATION ( cont inued ) 
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M(s) 

m b (s) 
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Effective mass moment of inertia of fluid 

Mass moment of inertia of disc 

Mass moment of inertia of mn slosh mass 

Mass moment of inertia of solidified fluid 

Imaginary unit 

Unit vector in x-direction 

Mass moment of inertia of fixed mass 

Effective mass moment of inertia of fluid 

Mass moment of inertia cf mn*'* 1 slosh mass 

s t 

Bessel function of the 1 kind of order v and argument 8 

First derivative of J (8) with respect to the 
argument 8 

Unit vector in y-direction 

Torsional stiffness coefficient of shaft 

Ratio of inner to outer radius of annular or annular- 
sector tank, b/a 

Unit vector in z-direction 

Pendulum length 

Fixed mass 

Magnitude of moment after n cycles 
Magnitude of initial moment 
Moment about q.-axis 

l 

Laplace transform of the moment M 
Laplace transform of the moment M on the baffle 
Mode of vibration in the ^-direction or y-direction 
Liquid mass 
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NOTATION (rout limed ) 


n s 1 osh mass 
Norma 1 direction 

Mode of vibration in the r -direct ion or x -direct ion 
Number of baffles 

Pressure in tbe q^-direction 
Gauge pressure 
Generalized coordinate 
Radius of circular cylindrical tank 
Cylindrical coordinate , (r ,<J» ,z ) 

Aspect ratio in xz-plane for rectangular tank, h/a 

Aspect ratio in yz-plane for rectangular tank, h/b 

Aspect ratio in xy-plane for rectangular tank, a/b 

Laplace transform variable 
Period of oscillation 
Time 

Potential function 

Timewise maximum fluid velocity 

A * 

Rotation vector for vorticity, £i + nj + 

Velocity component in x-direction 

AAA 

Velocity vector, ui + vj + wk 

Velocity component normal to tank wall 

Velocity of container 

Velocity component in y-direction 

Baffle width measured normal to tank wall 

Velocity component in z-direction 

Rectangular coordinate, (x,y,z) 

Displacement in x-direction 
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NOTATION (continued) 
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x(s) 

Y ( 6 ) 

V 

y'(p) 

V 


Translational displacement of n slosh mass 

Amplitude of displacement function, x = x e lajt or 

x = x sincjt 
o 

Laplace transformed rectangular coordinate 

Bessel function of the 2™* kind of order v and 
argument 3 

First derivative of Y v (3) with respect to argument (3 


y 


z 


Rectangular coordinate, (x,y,z) 
Translational displacement in y-direction 

Amplitude of displacement function, y = 

Rectangular coordinate, (x,y,z) 

Cyl indr i cal coordinate , (r , $ ,z ) 


Note: The notation given in articles 5.1.5.2 and 5. 1.7. 3 is applicable 

to these articles only. 
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Ratio of apex angle to 2 71 , &/ 2 V 

Baffle width parameter in Miles f equation, 2a /W 

Apex angle of sector tank or similar configuration 
Effective baffle area for a baffle located below 
the undisturbed free surface 

Effective baffle area for a baffle located above the 
undisturbed free surface 

Angular displacement of mn^ slosh mass 

Generalized coordinate, ^ ® 

n 

Roots of an eigenfunction where v = mn, n or (2m-l,n) 
Damping factor or damping ratio 

Surface displacement of liquid generally measured 
at tank wall 

Ratio of natural angular frequency to forcing frequency, 
w v /w, where v = mn, n, or (2m-l,n) 

Angular displacement of tank for pitching about y-axis 

iu>t 

Amplitude of displacement function, 0 = 0 Q e or 

6 = 0 sinwt 

o 

Laplace transform of the displacement 0 

Angular displacement of n^ pendulum 

Viscosity of liquid 

Mass density of liquid 

Total velocity potential, + $2 

Velocity potential of container 

Velocity potential of liquid 

Cylindrical coordinate, (r,<P,z) 

Angular displacement for roll about z-axis 

Amplitude of displacement function, 4> = or 

A = <j> sina)t 
T o 
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* 


Angular displacement for pitching about x-axis 
Amplitude of displacement function, X = X 
Angular displacement of disc relative to tank 


U) 


Forcing frequency 


r(v) 


Natural angular frequency for free liquid oscillations 
where v = mn, n, or (2m-l,n) 

Gamma function of argument v 
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SPECIAL DEFINITIONS 
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I. Introduction 


Fuel sloshing is defined as any periodic motion of a contained liquid 
propellant and results from the missile's oscillatory motion about its 
flight trajectory. The most likely causes of such disturbances are gust 
loads, control modes, and structural modes. If any of these excitations 

9 

have frequencies in the vicinity of the resonant frequency of a contained 
liquid, violent sloshing will occur. Since fuel and oxidizer amount to an 
extremely high percentage of the gross vehicle weight, the magnitudes of 
the liquid force and moment resultants are significant and cannot be 
neglected. A thorough knowledge of the magnitude and location of the 
external forces and moments acting on a space vehicle is required for 
making stability and structural investigations. Consequently, the effects 
, of liquid propellants on their containers while undergoing forced vibrations 
have been investigated both analytically and experimentally. The results of 

such studies are found in Chapter IV. 

The stability of a missile can be Increased by reducing the force and 

moment resultants caused by the oscillating propellant and Increasing the 
sloshing frequency. Increasing the eigen frequencies of the liquid de- 
creases the depth below the surface to which the disturbances extend. The 
above fact can be accomplished by utilizing longitudinal partitions which 
effectively reduce the sloshing liquid mass and Increase the frequencies of 
liquid oscillations. Also various types of baffles and floats can be used 
to help damp liquid motion. To evaluate the effectiveness of such anti- 
slosh devices, its damping factor and effect on the eigen frequencies of 
the fluid need to be known. Consequently, there has been a considerable 
amount of analytical and experimental work done in this area. The results 
of such investigations are found in Chapter V. 


1-1 . 


The exact formulation and solution to the problem of oscillations of 
a contained liquid with a free surface are extremely complex. To simplify 
the matter, it is convenient to assume that the fluid is nonviscous and 
the flow is irrotational , thereby permitting the use of potential theory. 
The problem then reduces to obtaining the velocity potential from Laplace's 
equation solved with the appropriate boundary conditions. 

The overall problem of missile stability can be investigated more 
easily if the liquid is replaced with a dynamically equivalent mechanical 
model. The equations of motion for a model are derivable from either La- 
grange's equations or simple equilibrium theory and from these equations, 
the force and moment resultants can be obtained. The characterist ics of 
the model are then determined so that the force and moment resultants, 
frequencies of oscillation, and mass and inertial characteristics are 
identically equal to those of the liquid. Also by using equivalent me- 
chanical models, the nonlinear nature of forced damped liquid oscillations 
can be approximated by equivalent linear damping. This makes it possible 
to obtain finite results near the resonant frequencies of the liquid, which 
was not possible in the fluid analysis, since the liquid was assumed to be 
inviscid. 
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IT. Fluid Theory (9) 


2.1 Assumptions 

The difficulty in formulating and obtaining the exact solution to the 
problem of oscillations of a contained liquid is simplified by making several 
assumptions concerning the nature of the fluid and the type of flow encount- 
ered. These assumptions are: (1) the fluid is both incompressible 
(p = constant) and frictionless (p = o), (2) the flow is irrotational 
(u = o), and (3) there are no sources or sinks; that is, the c6ntainer is 
neither draining nor filling. 

2.2 Basic Equations 

Irrotationality is a necessary and sufficient condition that the vel- 
ocity V can be expressed as the gradient of a velocity potential 4>, that is, 

= V (2-1) 


Because of incompressibility and the nonexistence of sources and sinks, the 
equation of continuity reduces to Laplace f s equation, 


V 2 # = U 


( 2 - 2 ) 


For an incompressible and frictionless fluid, irrotational flow, and con- 
servative body forces, those forces derivable from a potential function, 
Bernoulli's energy equation is 



+ U + = F(t) 


2-1 



where U is a potential function. Since gravity is the only important body 
force, then U = gz , and since the velocity is assumed to be small, the term 
containing V is neglected. Also since any arbitrary function of time can 
be added to the velocity potential without changing the flow it represents, 
the function F(t) can be conveniently absorbed by Consequently, the 
equation of motion reduces to 


r. + 82+ |i = ° 


From this expression, the pressure p is found to be 


(2-3) 


P = * p ( If + gZ \ 

For zero or uniform surface gauge pressure (p = 0 or p = p Q 
the surface displacements can be found as follows: 


(2-4) 
const • ) 


B| = {SS| .1 » 

surface g 3t 


(2-5) 
P, 


where, in the case of uniform surface gauge pressure, the constant — 

P 

is absorbed by The resultant liquid force acting on a container is then 
given by 


h L \ 


dA 


(2-6) 


A 


where is a generalized coordinate and and p^ are components of the 
force and pressure distributions in the q ^-direction. Taking the z-axis 
to be. the longitudinal axis of the container, the resultant liquid moment 
about a line perpendicular to the xz-plane and passing through the point 
(o,o,z) is given by 


H = f 

z=z J , 

(0,0, z) 


p (z + z)dA + 


x dA 


SIDE 


Bottom 


(2-7) 
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2 . 3 Boundary Conditions 


2.3.1 Velocity Potential. The total velocity potential ♦ is composed 
of two parts: 4^, the velocity potential of the container which is assumed 
to be small and 4^, the velocity potential of the liquid. These potentials 
must satisfy Laplace’s equation in the following manner: 

V 2 $ = V 2 (4> 1 + * ) = 0 (2-8) 

Since for a particular excitation the container velocity is known, the 
velocity potential of the container can easily be determined from equation 
(2-1) 

•■= Vj (2-9) 

where the resulting constant of integration can be absorbed by 4>^. 

2.3.2 Free Oscillations * Consideration of free oscillations is nec- 
essary in order to determine the natural frequencies of liquids in variously 
shaped containers. For the occurrence of free liquid oscillations, the 
container must be at rest (V^ = 0), thus 4>^ = constant. Consequently, 

4> = $2 since the constant can be absorbed by order to avoid separation 

over the wetted surface, the velocity of the fluid normal to the container 
wall must be identical to that of the container itself, that is, 


3*2 - 

— £ = v =0 

n 


( 2 - 10 ) 


Also, the pressure variation over the free surface must be equal to zero. 
Accordingly, the linearized free surface condition is 


3 *, 


( 2 - 11 ) 
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2.3.2 Forced Oscillations. When the container is subjected to harmonic 


excitation, the total velocity potential, ♦ = $ + , must be considered. 
Since over the wetted surface the normal velocity of the fluid must be the 
same as that of the container, it follows that 


3n 


■* 

V 

n 


(2-12) 


The linearized free surface condition corresponding to zero pressure var- 
iation at the free surface is 



a* 

9z 


0 


(2-13) 


The free surface boundary condition is actually composed of both a zero- 
pressure variation condition and a kinematic condition. In the kinematic 
condition, second and higher order terms were neglected thereby linearizing 
the free surface boundary condition. When the amplitude of the fluid is 
large relative to the tank dimensions, the effect of nonlinearity can become 
important, i*e., the magnitude of the higher order effects might not be 
negligible. Therefore, in order that the slope of the free surface be 
small, the resulting frequency of fluid oscillation must not be in the 
immediate neighborhood of a natural free surface frequency. Even so, 
relatively large amplitudes are possible when the container is also large. 
Also, it should be noted that the free surface condition given here and 
in Article 2.3.2 above is valid only in the case of a constant gravitational 
field. Consequently, the information presented in this handbook is applicable 
only for ground test applications. 

The expression for the natural frequencies for free liquid oscillations 
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is as follows: 



(2-14) 


^ e tanh (— e ) 
mn a mn a mn 

where 

g = the acceleration due to gravity 

a = the radius of the tank 

h = the depth of liquid in the tank 

e = the roots of an eigen function 
ran 

It is observed that the eigen frequencies are proportional to the square 
root of g. However, when the vehicle experiences a large acceleration 
normal to the mean free surface of the fluid, g becomes the effective 
longitudinal acceleration of the missile. Therefore, the natural frequencie.v 
of the liquid remain a constant while the vehicle is at rest but vary when 
the missile is in flight. It is also seen that the frequencies are inversely 
proportional to the square root of a and thereby decrease with an increase 
in the cross-sectional area of the container. The influence of tank geo- 
metry upon the frequencies of the liquid is also exhibited by the value of 
e mn » The liquid height has relatively no effect on the frequencies for 
h>a. As h decreases below a, the frequencies tend to zero. 

2.4 Initial Conditions 

Assume the harmonic forcing function is equal to zero for all time 
t equal to and less than zero (t < 0). 
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III. Equivalent Mechanical Model Theory 

An equivalent mechanical model is an assemblage of springs, dashpots, 
masses, and massless rods arranged in such a manner as to represent the 
dynamic behavior of a complex mechanical or non-mechanical system. Con- 
sideration is given here to those models which are dynamically equivalent 
to a liquid oscillating within its container. Dynamic equivalence is 
taken to mean the equivalence of force and moment resultants, frequencies 
of oscillation, and mass and inertial characteristics. 

There are two reasons for using equivalent mechanical models to repre- 
sent the sloshing behavior of a contained liquid. Since the problem of 
forced damped liquid oscillations is actually nonlinear in nature, an exact 
solution is practically impossible. However, through the use of an equiva- 
lent mechanical model, a good approximation can be made by introducing 
equivalent linear damping in the form of dashpots, thereby making it possible 
to obtain finite results near the resonant frequencies of the liquid. This 
was not possible in the fluid analysis, since the liquid was assumed to txe 
inviscid. Secondly, the overall problem of missile stability is simplified 
in that the vehicle’s equations of motion are not so complex. 

When a container partially filled with fluid is excited, the liquid 
in the bottom of the tank is little disturbed, whereas the liquid near the 
free surface oscillates. Analogously, the model is composed of a fixed mass 
in the bottom of the container to represent the essentially rigid or non- 
sloshing part of the liquid and a movable mass near the top of the container 
to represent the sloshing part of the liquid. The sum of the fixed and 
movable masses of the model are taken equal to the total mass of the liquid 
and the ratio of the movable mass to the total mass of the liquid increases 
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as the depth of the liquid in the container decreases. Also, the fixed mass 
and the movable mass (when the fluid is not oscillating) are located along 
the longitudinal axis of the container at such distances from the tank 
bottom as to yield the same inertial characteristics for rotation about any 
axis perpendicular to and passing through the longitudinal axis. Therefore, 
the fixed mass is located near the center of gravity of the undisturbed 
fluid and the location of the movable mass shifts towards the center of 
gravity of the fluid as the liquid level decreases. 

The kind of movable mass to be used and its restraints are dependent 
on the type of excitation encountered, namely: harmonic translation, 

harmonic pitching, or harmonic roll. For a particular model, harmonic 
translation and harmonic pitching have the following relationship: (a) 

translation refers to an excitation in a direction perpendicular to the 
longitudinal axis of the container while (b) pitching implies an excitation 
about an axis perpendicular to the direction of translation and the con- 
tainer^ longitudinal axis. For translational and pitching excitations, two 
possibilities exist to represent a movable mass and its constraints; either 
a spring-mass system whose motion is restrained by springs and dashpots 
aligned in the direction of translation or a pendulum-mass system having its 
massless lever arm attached to the longitudinal axis of the container. 
Because !f small oscillation theory 1 ’ has been assumed, the amplitude of the 
pendulum motion must be small, thereby maintaining the pendulum mass at a 
relatively constant distance from the container bottom. Therefore the 
inertial characteristics of the pendulum mass are essentially the same as 
those of the spring mass. 

For roll excitation, which is always taken to be about the longitudinal 
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axis of the container, two possibilities exist for representing s movable 
mass and its constraints; (a) a torsional pendulum whose shaft coincides 
with the container's longitudinal axis or (b) a torsional spring-mass system 
attached to the longitudinal axis of the container. 

Because of the defined arrangement of the constraints of movable 
masses with respect to the direction of excitation, the model becomes a 
single degree of freedom system, i.e., only one generalized coordinate 
is necessary to specify the position of the movable mass at any particular 
time. Note, however, that when a container partially filled with fluid is 
excited, the liquid can oscillate in various modes of vibration depending 
on the frequency and amplitude of the forcing function. Thus, for a 
circular cylindrical tank, or some other similar configuration, it is possible 
to have modes of vibration in both the radial and tangential direction. In 
other words, an oscillating liquid can exhibit several degrees of freedom. 

In order for a model to simulate these various modes of vibration, it is 
necessary to add additional movable masses to the model. Each new mass 
increases the model's degree of freedom by one and thereby makes it possible 
to represent an additional liquid mode of oscillation. For this reason the 
subscript mn is used on various parameters in order to clarify which mode of 
vibration is of interest, m is to represent the mth mode of vibration in the 
tangential direction and n is to represent the nth mode in the radial direct- 
ion for circular cylindrical tanks or similar configurations. In many cases, 
the only mth mode to be excited is nF=l. Thus, when only the subscript n is 
used, it is to be understood that nF=l . 

T « order for* » narticular movable mass to oscillate at the same 
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frequency as the liquid mode it represents, the following definitions have 
to be made: (1) the spring constant is defined as the slosh mass times 

the square of the natural angular frequency of the fluid, (2) the length 
of the pendulum axis is defined as the acceleration due to gravity over the 
square of the natural angular frequency of the fluid, and (3) the damping 
coefficient is defined as the slosh mass times the damping factor times the 
natural angular frequency of the fluid. Also, for the torsional pendulum, 
the torsional stiffness of the shaft is defined as the mass moment of inertia 
of the disc times the natural angular frequency and the damping coefficient 
of the disc is defined as the mass moment of inertia times the damping 
factor times the square of the natural angular frequency. The damping factor 
must be determined either analytically or experimentally (see Chapter V). 

Component values of the model are determined so that the forces and 
moments exerted by the model are identically equal to those exerted by the 
liquid. This is done by making a term-wise comparison of the force and 
moment equations derived from hydrodynamic aspects with those derived from 
the model. In order to put the fluid and model equations into similar form 
for comparison, it is many times necessary to expand certain terms of these 
equations into series. Even though the model equations contain damping 
terms and the fluid equations do not, the comparison of equations is still 
valid because the modes of vibration change little as long as the amount 
of damping remains small. 
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IV. Rigid Tanks 


4.1 Introduction 

The following sections present results of analytical studies of liquid 
oscillations in containers of various shapes when subjected to different 
types of excitation. Each section is concerned with a particular container 
shape, either circular cylindrical, sector, quarter-sectored, eighth- 
sectored, annular, annular-sector, or rectangular and is preceded by a 
table of contents, a diagram showing the container and its coordinate system, 
and other information pertinent to the material that follows* The containers 
are assumed to have exhibited rigid-body behavior* Th e boundary conditions 
and resulting fluid velocity potentials, natural frequencies and liquid 
force and moment resultants are given along with equivalent mechanical 
models. The elements of these models are defined both analytically and 
graphically. 

The orthogonal coordinate systems to which the force and moment result- 
ants are referenced are identical to those shown in the figures, the coordin- 
ate origins always being located at the center of gravity of the undisturbed 
fluid. However, in the fluid analysis, it is frequently convenient to 
translate the coordinate origin along the z-axis to the free surface of the 
undisturbed liquid, thereby referencing the velocity potential to such a 
system. This fact presents no difficulty since equation (2- 7 ), used for 
calculating the resultant moment, gives the moment referenced to any point 
(o,o,z). Therefore, it is necessary to take z = -h/2 when the velocity 

potential is referenced to the fluid surface and z = 0 when the velocity 

potential is referenced to the center of gravity of the undisturbed fluid. 


The orthogonal coordinate system to which a particular velocity 



potential is referenced is determined in the following manner. The z-axis 
is always identical to that shown in the appropriate figure. The location 
of the coordinate origin along the z-axis will be either at the center of 
gravity of the undisturbed fluid or at the undisturbed free liquid surface. 
Upon examination of the z -coordinate at which the corresponding free surface 
boundary condition is evaluated, the location of the coordinate origin is 
obvious . 

Also, it should be mentioned that, for a particular container, super- 
position of force and moment resultants resulting from several types of 
excitations is possible because of the linearized theory used herein. 
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4.2 Circular Cylindrical Tank 
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Container: The tank is a right circular cylinder of radius a and is filled 

with a liquid to a depth h. 

Coordinate System : The origin of the xyz-system is located at the center of 
gravity of the undisturbed fluid and the origin of the x^x 1 system is 
located at the geometric center of the tank base# 

References : (4, 44, and 52) 

Comments : The terms I and c, occurring in the model elements, must be 

experimentally determined, generally from torsion spring experiments 
of a sealed container of fluid. Their values are dependent upon 
container shape and the type of liquid. 
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Table 4-1. Boundary Conditions, Velocity Potential, and Natural Frequency (continued) 


o 
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where e are roots of J,(e ) = 0. (see Table A-3) 


Circular Cyli ndrical Tank | Excitation: Harmonic Translation or Pitc hing 

Translation in the x-direction, x = x e ltut . 
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Table 4-3. Model Analysis 


Circular 


Cylindrical Tank 


Spring-Mass Model 


Excitation : Harmonic Translation and/or Pitching 


Figure 4-1 shows a diagram of the spring-mass model used in representing 
the dynamic response of a liquid in a circular cylindrical tank when sub- 
jected to harmonic translation in the x-direction and/or pitching about the 
y-axis . 

Coordinate System: 

The origin is located at the center of gravity of the undisturbed 
liquid. 


Model Description : 

The components of the system are as follows: 

1. A fixed mass M having a moment of inertia I is rigidly connected 
to the tank and is located on the z-axis at a distance H below 
the coordinate origin. 

2. A set of movable masses m is distributed along the z-axis when 
the tank is at rest at distances h above the origin. These modal 
masses are constrained by spring-dashpot systems, having spring 
stiffness coefficients k and viscous damping coefficients c , 

to remain in the xy-plane and to move only in a direction parallel 
to the x-axis. Translational displacements of these masses with 
respect to the container are denoted by x^. 

3. A massless disc having a moment of inertia I, is located at the 
coordinate origin. Its motion is confined to rotation about the 
y-axis and is cushioned by a dashpot having a viscous damping 
coefficient c,. The angular displacement of the disc relative to 
the tank is defined by 

Equations of Motion: 

The equations, obtained through Lagrange 1 s equations, are as follows: 


1. Force Equation: 


F 

x 


M(x + HO) 


l m (x + x + h e) 
n n n° 


2. Moment Equation: 


M = - 

y 


- I d (e + $) + g 


00 


l 

n=l 


m x 
n n 


00 


I 

n=l 


m h 
n n 



+ h n ° 5 


k-Q 


(1 + MH 2 )6 






Table 4.3 Model Analysis (continued) 


Circular Cylindrical Tank 


Spring-Mass Model 


Excitation: Harmonic Translation and/or Pitching 


Equations of Motion (continued): 

3. Disc Equation: 

I d (9 + *) + cj> = 0 

4. Slosh -Mass Equation: 

m (x + x +h ®) + mg ux +kx - m g0 = 0 
n n n n°n n n n n n° 

From these equations, the model force in the x-direction and the moment about 
the y-axis can be found (see Table 4-4). 
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Circular Cylindrical Tank | Spring-Mass Model 

Excitation; Harmonic Translation and/or Pitching 
Figure 4-1. Equivalent Mechanical Model 
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Table 4-5. Model Elements 


Circular Cylindrical Tank J Spring -Mass Model 


Excitation: Harmonic Translation and/or Pitching 


Natural Frequency u 2 & e tanh k 

nan n 


(fig. 4-2, 4-3) 


Spring Constant 


k in u) 
n n n 


(fig. 4-4, 4-5) 


Damping Coefficient 
of Slosh Mass 


Damping Coefficent 
of Disc 


c m g w 
n n n n 



(fig. 4-6, 4-7) 


Ratio of Slosh Mass 
to Fluid Mass 


in 2 tanh (< ) 
n n__ 

L ( c - 1 ) < 
n n 


(fig. 4-8) 


Ratio of Fixed Mass M mm 

to Fluid Mass — 1 - T -p- 

m, L m_ 


(fig. 4-9, 4-10) 


Ratio of Slosh Mass | |^ n | ^ 

Coordinate to Fluid 
Depth 


Ratio of Fixed Mass 
Coordinate to Fluid 
Depth 


Moment of Inertia 
of Solidified Fluid 


Moment of Inertia 
of Disc (c^ ~ 0) 


Moment of Inertia 
of Disc (c^ f 0) 


Moment of Inertia 
of Fixed Mass 



(fig. 4-12) 


2 • [ 1 " (2/k ) tanh (x /2)] 

r L n n J 


I, / v 2 • 1 - (2/te ) 

= s }) I 1 r ^- 

n T h 2 V7 n=o (c 2 - 1) 



(fig. 4-13) 
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Ratio of Fluid Depth to Tank Radius 
Figure 4-3 Model Element Graph 
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FLUID HEIGHT RATIO 
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CIRCULAR CROSS SECTION 
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F-fonre 4-5. Model Element Graph 
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Ratio of Fluid Depth to Tank Radius 
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Figure 4-9. Model Element Graph 
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Figure 4-10. Model Element Graph 
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Figure 4-11. Model Element Graph 



MOMENT OF INERTIA OF THE DISC 
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Figure 4-12. Model Element Graph 
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Table 4-6. Model Analysis 


Circular Cylindrical Tank 


Pendulum Model 


Excitation: Harmonic Translation and/or Pitching 


Figure 4-14shows a diagram of the pendulum model used in representing the 
dynamic response of a liquid in a circular cylindrical tank when subjected 
to harmonic translation in the x-direction and/or pitching about the y-axis. 


Coordinate System: 

The origin is located at the center of gravity of the undisturbed liquid. 


Model Description : 

The components of the system are as follows: 

1. A fixed mass M having a moment of inertia I is rigidly connected to 
the tank and is located on the z-axis at a distance H below the 
coordinate origin. 

2. A set of movable masses m is distributed along the z-axis when the 
tank is at rest. These modal masses are pendulums having massless 
lever arms of length L attached to the z-axis at distances H above 
the origin. They are constrained by dashpots having viscous 3amping 
coefficients c n to remain approximately in the xy-plane and to move 
approximately parallel to the x-axis. Angular displacements of the 
pendulum with respect to the tank (z-axis) are denoted by X^. 

3. A massless disc having a moment of inertia 1^ is located at the 
coordinate origin. Its motion is confined to rotation about the 
y-axis and is cushioned by a dashpot having a viscous damping 
coefficient c^. The angular displacement of the disc relative to 
the tank is defined by ip. 


Equations of Motion: 

The equations, obtained through Lagrange 1 s equations, are as follows: 
1. Force Equation: 

F = - M(x + H0) - l m f x + L X + (H - L )8 1 
x =i n L nn nn j 


2. Moment Equation: 

OO 00 

M = - (I + MH 2 )e - 1,(0 + $) + g Y tn L A - Y m(H - L ) • 
y d ° L ,nnn‘‘,nnn 


n=l 


n=l 


LA + (H - L )0 
n n n 
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Figure 4-6. Model Analysis (continued) 


Circular Cylindrical Tank 


Pendulum Model 


Excitation: Harmonic Translation and/or Pitching 


Equations of Motion (continued): 


3. Disc Equation: 


1^ (6 + + c^ = 0 

4. Slosh-Mass Equation: 


m x+LX + (H -L)6 
n n n n n 


mgwLX - mg 
nn n n n n 


6 = 0 


From these equations, the model force in the x-direction and the moment about 
the y-axis can be found (see Table 4-7). 
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Circular Cylindrical Tank 


Pendulum Model 


Excitation: Harmonic Translation and/or Pitching 

Figure 4-14. Equivalent Mechanical Model 
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Table 4-8. Model Elements 


Circular Cylindrical Tank 


Pendulum Model 


Excitation: Harmonic Translation and/or Pitching 


2 o 

Natural Frequency w = e tanh < 


nan 


(fig. 4-15) 


Pendulum Length 


L = g/< 

n n 


Damping Coefficient 
of Slosh Mass 


c = m g a) 
n n n n 
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c = c 1 + * y 

d l » 2 <i - i ) 2 

s 


„ - 01 , w m 2 tanh (k ) 
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"L (e 2 - 1) tc 
n n 
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= 1 - I - 

l* TTI 


"L 


Ratio of Pendulum H - L 
4/ ^ I 1 n n 

Mass Coordinate 
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of Solidified Fluid 


(fig. 4-16, 4-17) 
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SUSPENSION POINT LOCATION 
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Figure 4-15. Model Element Graph 
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Table 4-9. Model Analysis 


Circular Cylindrical Tank 


Pendulum Model 


Excitation: Arbitrary Translation and/or Pitching 


Figure 4-18 shows a diagram of the pendulum model used in representing the 
dynamic response of a liquid in a circular cylindrical tank when subjected to 
an arbitrary translation in the x-direction and/or pitching about the y-axis. 

Coordinate System: 

The origin is located at the geometric center of the tank bottom. 


Model Description : 

The components of the system are as follows: 

1. A fixed mass M having a moment of inertia I is rigidly connected to 
the tank and is located on the z-axis at a distance H above the 
coordinate origin. 

2. A set of movable masses m n is distributed along the z-axis when 

the tank is at rest. The modal masses are pendulums having massless 
lever arms of length 1 ^ attached to the z-axis at distances H n above 
the origin. Angular displacements of the pendulums with respect to 
the tank (z-axis) is denoted by X n and 0 is a space fixed coordinate 
defining the pitching angle of the container. Angular displacements 
of the pendulums are therefore defined by a generalized space fixed 
coordinate, r n = X n " 


Equations of Motion: 

The equations, obtained through Lagrange f s equations, are as follows: 

1. Force Equation: 

F = -(M-hm)x-lMH+m(H - L )1 0 - m LX 

n L nn nJ nnn 

2. Moment Equation: 

M = - f MH + m (H - L )1 x 

L n n n J 

+ [ MH + m (H - L )] 0 
L n n n J 

+ m gL X 
n° n n 

3. Slosh Mass Equation: 

X + o)^X = - 7 — I x + (H - L )0 - g0 
n n \ ^ nn 6 J 

Using the method of Laplace transforms, the force and moment have been repre- 
sented as transformed variables, which thereby become a function of the La- 
place-transform operator s (see Table 4-10). 


- I MH 2 + m (H - L ) 2 + l( 6 

L nnn J 

- m L (H - L )A 
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Excitation: Arbitrary Translation and/or Pitching 


Figure 4-18. Equivalent Mechanical Model 
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Circular Cylindrical Tank Pendulum Model 
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Table 4-11. Model Elements 


Circular Cylindrical Tank Pendulum Model 

Excitation: Arbitrary Translation and/or Pitching 

Nat iona 1 Frequency 

2 o 

a) = ■* e tanh k 
nan n 

Pendulum Length 

L =& 2 
n 

n 

Slosh Mass 

w 

2 tanh k 

n 

.< ( e 2 - 1). 
n n 
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00 

M = “l ■ E m n 
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|h - L 1 = h 
n n 

2 - cosh < 
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n n J 
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I I L 1 a n n n 

11 m T -m 2 /t 2 2,2 . . . 

LnL 4h tc(e - U cosh k 

n n n 
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i 

3 

2 “ 
8k +(2k -8)sinhK -2 k coshK 
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4.3 Sector Tank 
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Container : The tank is a sector of a right circular cylinder of radius a 

and filled with a liquid to a depth h. The vertex angle is denoted 
by a. 

Coordinate System : The origin is located at the center of gravity of the 

undisturbed fluid which would be contained in a right circular cylinder 
generated by revolving the sector tank about the z-axis. The x-axis 
must lie in the sector wall. 


References : (4) 

IT 3 TT 

Comments : The results given in this section are not valid for a = — — 
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Table 4-12* Boundary Conditions, Velocity Potential, and Natural Frequency 


o 
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where e are roots of J /0 (e „ 
mn m/2a mn 



Table 4-12. Boundary Conditions, Velocity Potential, and Natural Frequency (continued) 
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where e are roots of J , ' 

mn m/2a 
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Table 4-14. Boundary Conditions, Velocity Potential, and Natural Frequency 
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Natural angular frequency: 


lahle 4-14. Boundary Conditions, Velocity Potential, and Natural Frequency (continued) 



where e are zeros of J /0 (e ^ 
mn m/2a mn 
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Table 4-15. Liquid Force and Moment Resultants (continued) 
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Table 4-16. Boundary Conditions, Velocity Potential, and Natural Frequency 




Table 4-17. Liquid Force and Moment Resultants 
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Table 4-17. Liquid Force and Moment Resultants (continued) 
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.■ontalner ; The tank is a right circular cylinder of radius a divided by 

longitudinal partitions intersecting at right angles, a — it/ 2, and is 
filled with liquid to a depth h. 

Coordinate System : The origin is located at the center of gravity of the 

undisturbed fluid and the x-axis must lie in a sector wall. 

References : (4, 60) 

Comments : The results given in this section are applicable to only one 

sector, i.e., the boundary conditions and resulting velocity potential, 
natural frequency, and force and moment resultants are responses due to 
liquid oscillations in only one sector (except for the material concern- 
equivalent mechanical models which applies to the complete tank of 
four sectors). 
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(see Table A-4). 


Quarter-Sectored Tank | Excitation: Harmonic Pitchin 
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are roots of J 0 (e ) — 0. (see Table A-4) 



Table 4-19. Liquid Force and Moment Resultants 
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Table 4-20. Model Analysis 


Quarter-Sectored Tank | 

Excitation: Harmonic Translation 


Spring-Mass Model 
and/or Pitching 


Figure 4-19 shows a diagram of the spring-mass model used in representing 
the dynamic response of a liquid in a quarter-sectored tank when subjected 
to harmonic translation in the x-direction and/or pitching about the y-axis. 


Coordinate System: 

The origin is located at the center of gravity of the undisturbed liquid. 


Model Description: 

The components of the system are as follows: 

1. A fixed mass M having a moment of inertia I is rigidly connected 
to the tank and is located on the z-axis at a distance H below 
the coordinate origin. 

2. A set of movable masses m is distributed along the z-axis when 
the tank is at rest at distances h above the origin. These modal 
masses are constrained by spr ing-dashpot systems, having spring 
stiffness coefficients k and viscous damping coefficients c , 

to remain in the xy-plane n and to move only in a direction parallel 

to the x-axis. Translational displacements of these masses with 

respect to the container are denoted by x 

J mn 

3. A massless disc having a moment of inertia 1^ is located at the 
coordinate origin. Its motion is confined to rotation about the 
y-axis and is cushioned by a dashpot having a viscous damping 
coefficient c^ . The angular displacement of the disc relative to 
the tank is defined by ip. 


Equations of Motion: 

The equations, obtained through Lagrange T s equations, are as follows: 
1. Force Equation: 

OO 00 

F = - M(x + H0) - y y m (x + x + h 0) 
x u u - mn mr mn 

m=l n=T 


2. 


Moment 
M = - 

y 


Equation : 

OO 00 

(I + MH 2 )e - I H (0 + 5) + g y y m x 

d 6 L n L t mn mn 

oo oo m=1 n=1 

“I 7 m h (x +h0) 

, u - mn mn mn mn 
m _ ~l n=l 
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Table 4-20. Model Analysis (continued) 


Quarter -Sectored 


Tank 


Spring-Mass Model 


Excitation: Harmonic Translation and/or Pitching 


Equations of Motion (continued): 

3. Disc Equation: 

I d (5 + ♦) + c d * = 0 

4. Slosh-Mass Equation: 

m (x + x + h6)+mgux +kx -mg0-O 
mn mn mn mn mn mn mn mn mn mn 

From these equations, the model force in the x-direction and the moment about 
the y-axis can be found (see Table 4-21). 
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Figure 4-19. Equivalent Mechanical Model 
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Table 4-22. Model Elements 


Quartered-Sectored Tank 


Spring-Mass Model 


Excitation: Harmonic Translation or Pitching 


Natural Frequency 


^ e tanh < 
mn a mn mn 


Spring Constant 


k = m a) 
mn mn mn 


Damping Coefficient 
of Slosh Mass 


c — m g a) 
mn mn mn mn 


64 tanh tc 


mn 


Ratio of 
Slosh Mass 
to Fluid Mass 


m 

mn = 

m T 2 ,2 .2x2 

L n < U -4m )J 0 

mn mn 2m mn 


+ 2 l J 2m+2U+I (e mn ) 


U=o 


1 00 J 

l — 

J u=o ( 2m 


G mn^2n/ e mn^ (fl-S* ^“20, 

4-21) 


(4m -1) 


2m+2u+l (g mn ) 


(2m+2u+3)(2m+2u-l) 


Ratio of Fixed 
Mass to Fluid Mass 


oo « m 

r r mn 

m. ^ ^ til. 

L ri=o m=o L 


i 


Ratio of Slosh 
Mass Coordinate 
to Fluid Depth 


mn I 
h 


tanh 


mn 


mn 

2 


(fig. 4-22) 


Ratio of Fixed 
Mass Coordinate 
to Fluid Depth 


ni T 


00 / m 


lUl = _L y y 

h M L L 

m=o n—o 


mn ' 


/*. 


m T 


mn 
h j 


Moment of Inertia 
of Solidified 
Fluid 


m L h‘ 


= (^-) 2 [— (ll 
n ' 12 V 


)2 + i 


Moment of Inertia 


m L h 


= ( 


0 f . * “ ( - 1 ) m a b 

f) 2 U+8 l l ^ 


m^o n==o 7T a c 


mn 


i 2 ^ , mn 

1 - — tanh — 

mn 




















SLOSHING MASS 
VS 

FLUID HEIGHT RATIO 
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Figure 4-20. Model Element Graph 






Figure 4-21. Model Element Graph 
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MASS LOCATION HEIGHT RATIO 

VS. 

FLUID HEIGHT RATIO 



o 2 4 e e to 

Ratio of Fluid De F :h to T* n k Radius 
Figure 4-22. Model Element Graph 
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Table 4-23. Model Analysis 


Quarter-Sectored Tank 


Pendulum Model 


Excitation: Harmonic Translation and/or Pitching 


Figure 4-23 shows a diagram of the pendulum model used in representing the 
dynamic response of a liquid in a quarter-sectored tank when subjected to 
harmonic translation in the x-direction and/or pitching about the y-axis. 


Coordinate System: 

The origin is located at the center of gravity of the undisturbed liquid. 


Model Description : 

The components of the system are as follows: 

1. A fixed mass M having a moment of inertia I is rigidly connected to 
the tank and is located on the z-axis at a distance H below the 
coordinate origin. 

2. A set of movable masses m is distributed along the z-axis when the 
tank is at rest. These mo3al masses are pendulums having massless 
lever arms of length L attached to the z-axis at distances H above 
the origin. They are constrained by dashpots having viscous damping 
coefficients c to remain approximately in the xy-plane and to move 
approximately parallel to the x-axis. Angular displacements of the 
pendulum with respect to the tank (z-axis) are denoted by X^ . 

3. A massless disc having a moment of inertia 1^ is located at the 
coordinate origin. Its motion is confined to rotation about the 
y-axis and is cushioned by a dashpot having a viscous damping 
coefficient c^. The angular displacement of the disc relative to 
the tank is defined by ip. 


Equations of Motion: 

The equations, obtained through Lagrange's equations, are as follows: 
1. Force Equation: 


F 

x 


- m(x + lie) 


CT> CO — 

y y m x + l 

is ran l ran 

m— I n— 1 



L 

mn 



2. Moment Equation: 


M ~ - (I H“ MH 2 )6 - 1.(0 + *) +g y y m L X 

V d T ° ** L . mn mn 


m=l n=l 


mn mn mn 


l l 

m— 1 n— 1 


m (H - L ) L X +(H-L)0 
mn mn mn mn mn mn 
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Table 4-23. Model Analysis (continued) 


Quarter -Sectored Tank 


Pendulum Model 


Excitation: Harmonic Translation and/or Pitching 


Equations of Motion (continued): 

3. Disc Equation: 

1^ (0 + ifO + = 0 

4. Slosh-Mass Equation: 


m 


mn 


x + L X + (H - L )0 
mn mn mn mn 


+ m g w L A -mg 
mn mn mn mn mn mn 


e = o 


From these equations, the model force in the x-direction and the moment about 
the y-axis can be found (see Table 4-24). 








Quarter-Sectored Tank 


Pendulum Model 


Excitation: Harmonic Translation and/or Pitching 


Figure 4-23. Equivalent Mechanical Model 
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Translation in the x-direction 
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Table 4-25. Model Elements 


Quarter-Sectored Tank I Pendulum Model 


Excitation: Harmonic Translation and/or Pitching 


= £ r t-arVh yc 


Natural Frequency w = * c tanh < 

M J mn a mn mn 


Pendulum Length 



Damping Coefficient 
of Slosh Mass 


c = m g in 
mn ran mn mn 


Ratio of 
Slosh Mass 
to Fluid Mass 


64 tanh < 


A (e 2 -4m 2 )j 2 (e ) L (4m 2 -l) 
mn mn 2m mn L 


E *1 o ( £ ) 

mn zm mn 


+ 2 ^ J 2«+ 2 UH (e n,„ ) £ 

y=o _ y=o 


(e ) 


2m+2u+l mn 
(2m+2u+3)(2m+2u-l) 


Ratio of Fixed 
Mass to Fluid Mass 


Ratio of Pendulum 
Mass Coordinate to 
Fluid Depth 


Ratio of Fixed 
Mass Coordinate 
to Fluid Depth 


Moment of Inertia 
of Solidified 
Fluid 


Moment of Inertia 

of Disc (c, = o) 
a 


M_ _ y y mr 

m. ^ ^ m, 

X n=o L 



- <s> 2 fc + i 


I d ( a \2 { 1 . n r a m b nm If. 2 t . 'jnn 

— 2 = <h } 2 + 8 l 2 {'-T-*™*— 

m T h m — o n — o xac \ mn 

L V. mn w 


J 9tn^ e mn ^ \ 2C 2 L (e ) tanh < 

2m mn , „ „ \ mn 2 mn mn 

2 + (e >J - ; 

° ™ (fa n 2 -l) k 

1 mn 
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Table 4-26. Boundary Conditions, Velocity Potential, and Natural Frequency 
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Table 4-27. Liquid Force and Moment Resultants 
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Table 4-27. Liquid Force and Moment Resultants (continued) 
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Table 4-28* Model Analysis 


Quarter -Sectored 


Tank 


Torsional Pendulum Model 


Excitation: Harmonic Roll 

Figure 4-24 shows a diagram of the torsional pendulum model used in 
representing the dynamic responseof a liquid in a quarter-sectored tank when 
subjected to harmonic roll about the z-axis. 


Coordinate System: 

The origin is located at the center of gravity of the undisturbed liquid- 


Model Description: 

1. A fixed mass having a mass moment of inertia J is rigidly connected 
to the tank and is located on the z-axis at the base of the tank. 

2. A disc having a mass moment of inertia J mn is connected to a shaft 
having a torsional stiffness K mn which is in turn rigidly attached 
to the tank. Viscous damping is introduced along the periphery 

of the disc, the damping coefficient of which is c mn * Angular 
displacement of the disc relative to the tank is denoted by 6 mn * 


Equations of Motion: 

The equations, obtained i“hrough Lagrange r s equations, are as follows: 
1. Moment Equation: 


« " - J 'in + I 


m—1,3,5,... n=o 
2. Slosh Mass Equation: 


l J (♦ + B ) 

u mn 


mn 


l 


m— 1 ,3,5,..* n-o 


l j (♦ + 8 ) + y 

L mn mn L 


m=l ,3,5,... n=o 


l c 6 

mn i 


mn mr 


I l K 6 = 0 

0 c _ mn mn 

m n—o 


From these equations, the model moment about the z-axis can be found 
(see Table 4-29). 
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Quarter -Sectored Tank 


Torsional Pendulum Model 


I 

Excitation: Harmonic Roll 

Figure 4-24- Equivalent Mechanical Model 
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Table 4-29. Moment Resultants 











Table 4-30. Model Elements 


Quartered -Sectored Tank Torsional Pendulum Model 

Excitation: Harmonic Roll 


Natural Frequency 

w 2 = f e tanh k 

mn 2 mn mn 

Torsional Stiffness 
of Shaft 

K = J to 2 
mn mn mn 

Damping Coefficient 
of Slosh Disc 

II 

r 

Moment of Inertia 

of Slosh Disc 

(c = o) 
mn 

8m. a 2 T (2f -e ) L (e ) tanh * 

_ L ( n n I n n t 

"' m » 2 T K n 

2( V”-> L l <C m > tanh Vnl 


m(m^-l) < J 

mn 

Moment of Inertia 

of Slosh Disc 

(c f o) 

mn 

8m.a 2 (n 2 +ig n -1) f (2f -e ) L (e ) tanh < 

L mn ®mn mn 1 n n 1 n n 

mn *2 [ 

2 ( f -me ) L (e ) tanh * ] 

mn mn 1 mn mn 

(n n'° K n 
n n 

m 2 (m 2 -l)(ri 2 -1) < 

mn mn 

Moment of Inertia 
of Solidified Fluid 

T - 1 2 

J s 2 m L a 

Effective Moment of 
Inertia of Fluid 

j = J fl - % l -r 1 J 

s L ir nf=l ,3,5... m (m+1) j 


Moment of Inertia of 
Fixed Mass 

00 00 

j = 3 - y I j 

m— 1 ,3,5... n~o 












Slosh Mass Moment of Inertia 
vs Fluid Height Ratio for a Quarter 
Sector Tank When Subjected 
to Harmonic Roll 
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Figure 4-25. Model Element Graph 



Ratio of Slosh Mass Moment of Inertia 
to Fluid Mass Moment of Inertia for 
a Quarter Sector Tank When 
Subjected to Harmonic Roll 


co 
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Figure 4-26. Model Element Graph 



4.5 Eighth-Sectored Tank 
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Container : The tank is a right circular cylinder of radius a divided by 

longitudinal partitions intersecting at equal angles, a = tt/ 4, and is 
filled with liquid to a depth h. 

Coordinate System ; The origin is located at the center of gravity of the 
undisturbed fluid and the x-axis must lie in a sector wall* 

References t (60) 

Comments : The results given in this section are applicable to only one 

sector, i.e., the boundary conditions and resulting velocity potential, 
natural frequency, and force and moment results are responses due to 
liquid oscillations in only one sector (except for the material concern- 
ing equivalent mechanical models which applies to the complete tank of 
eight sectors). 
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Boundary Conditions, Velocity Potential, and Natural Frequency 
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Natural angular frequency: 


Liquid Moment Resultant 











Table 4-33. Model Analysis 


Eighth-Sectored Tank 


Torsional Pendulum Model 


Excitation: Harmonic Roll 


Figure 4-7 shows a diagram of the torsional pendulum model used in 
representing the dynamic response of a liquid in an eighth-sectored tank when 
subjected to harmonic roll about the z-axis. 


Coordinate System: 

The origin is located at the center of gravity of the undisturbed liquid. 


Model Description: 

1. A fixed mass having a mass moment of inertia J is rigidly connected 
to the tank and is located on the z-axis at the base of the tank. 

2. A disc having a mass moment of inertia is connected to a shaft 

having a torsional stiffness K^ n which is in turn rigidly attached 
to the tank. Viscous damping is introduced along the periphery 

of the disc, the damping coefficient of which is c rnn . Angular 
displacement of the disc relative to the tank is denoted by 


Equations of Motion: 

The equations, obtained through Lagrange r s equation, are as follows: 


1. 

Moment Equation: 





00 

CO 



M = 
z 

- JB +1 

mn o c 

m i,_j,3,..« 

l 

n=r> 

J (♦ + 6 ) 

mn mn 

2. 

Slosh 

Mass Equation: 




oo 

00 


00 00 


I 

mr=l,3 

[ j (♦ + 

_ * mn r 

,0 , . . . n—o 

6 ) 
mn 

+ y y c 8 

L , c L mn mn 

m=l ,3,5,... n^ 


00 


+ l 

m— 1 ,3,5 


1 

n=o 


K 6 

mn mn 


- 0 


From these equations, the model moment about the z-axis can be found 
(see Table 4-34). 
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Eighth-Sectored Tank 


Torsional Pendulum Model 


Excitation: Harmonic Roll 


Figure 4-27. Equivalent Mechanical Model 
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Table 4-34* Model Moment Resultant 












Table 4-35. Model Elements 


Eighth-Sectored Tank Torsional Pendulum Model 

Excitation: Harmonic Roll 

Natural Frequency 

2 e 

co = & e tanh k 

mn a mn mn 

Torsional Stiffness 
of Shaft 

K = J u> 2 

mn mn mn 

Damping Coefficient 
of Slosh Disc 

c = J g w 

mn mn mn mn 

Moment of Inertia 

of Slosh Disc 

(c = 0 ) 

mn 

32m a 2 f -(f -2me ) L (e ) tanh k 1 (fig. 
L | mn mn 2 mn mn | v . 

J n 1 <-) 4"28j 

run 71 L m(4m -1) tc J 4-29) 

mn 

Moment of Inertia 

of Slosh Disc 

(c f 0 ) 

mn 

-32m a ^ (n ^ + ig n - 1) 

L mn mn mn 

mn 2, 2 . * 

it (n - 1) 

mn 

(f -2me ) L (e ) tanh < 
mn mn 2 mn mn 

m(4m^ - 1) tc 

mn 

Moment of Inertia 
of Solidified Fluid 

J s = i v 2 

Effective Moment of 
Inertia of Fluid 

5 “ J s l 1 ‘ 2 ^ - 2, ,2 ] 

L it m=l ,3 ,5 , . • • m (m+1) J 

Moment of Inertia 
of Fixed Mass 

00 00 

j = j - 1 1 j 

0 c ^ mn 

j j . . . n — 0 
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Slosh Mass Moment of Inertia 
vs Fluid Height Ratio for an Eighth 
Sector Tank When Subjected to 
Harmonic Roll 
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Figure 4-28. Model Element Graph 








Ratio of Slosh Mass Moment of Inertia 
to Fluid Mass Moment of Inertia for 
an Eighth Sector Tank When Subjected 
to Harmonic Roll 



Figure 4-29. Model Element Graph 
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Container : The tank consists of two concentric right circular cylinders 

having an inner radius r = b and an outer radius r = a. The space be- 
tween the concentric walls is filled with liquid to a depth h. 

Coordinate System : The origin is located at the center of gravity of 
the undisturbed fluid. 

References : (52 and 63) 

Comments : The terms I and c, occurring in the model elements, must be 

experimentally determined, generally from torsion spring experiments 
of a sealed container of fluid. Their values are dependent upon 
container shape and the type of liquid. The model element graphs 
are given for k = b/a = 0.5. These graphs are, within +10 per cent, 
applicable for values 0.3 < k < 0.7. 
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Table 4-36. Boundary Conditions, Velocity Potential, and Natural Frequency 
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Table 4-36. Boundary Conditions, Velocity Potential, and Natural Frequency (continued) 
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are roots of A.U ) = J,(£ )Y.'(k£ ) - J.'(k£ ) Y '<£ ) = 0. (see Table A-6) 
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Table 4-37. Liquid Force and Moment Resultants 
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Table 4-38. Model Analysis 


Annular Tank 


Spring-Mass Model 


Excitation: Harmonic Translation and/or Pitch ing 

Figure 4-30shows a diagram of the spring-mass model used in representing 
the dynamic response of a liquid in an annular tank when subjected to harmonic 
translation in the x-direction and/or pitching about the y-axis. 


Coordinate System: 

The origin is located at the center of gravity of the undisturbed 
liquid . 


Model Description: 

The components of the system are as follows: 

1. A fixed mass M having a moment of inertia I is rigidly connected 
to the tank and is located on the z-axis at the distance H below 
the coordinate origin. 

2. A set of movable masses m^ is distributed along the z-axis when 
the tank is at rest at distances h above the origin. These modal 
masses are constrained by spring-dashpot systems, having spring 
stiffness coefficients k and viscous damping coefficients c , 

to remain in the xy-plane and to move only in a direction parallel 
to the x-axis. Translational displacements of these masses with 
respect to the container are denoted by x^. 

3. A massless disc having a moment of inertia 1^ is located at the 
coordinate origin. Its motion is confined to rotation about the 
y-axis and is cushioned by a dashpot having a viscous damping 
coefficient c^. The angular displacement of the disc relative to 
the tank is defined by 


Equations of Motion: 

The equations, obtained through Lagrange’s equations, are as follows: 

1. Force Equation: 

oo 

F = - M(x + H0) - V m (x + x + h *0) 
x L , n n n 

n=l 

2. Moment Equation: 

co 

m = - (i + mh 2 )0 - i J (e + ip) + g y 

y d n=l 
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Table 4-38. Model Analysis (continued) 


Annular Tank 


Spring -Mass Mode 


1 


Excitation: Harmonic Translation and/or Pitching 


Equations of Motion (continued): 
3 Disc Equation: 

T (6 + l) + c J = o 0 

d d 

4. Slosh-Mass Equation: 


m (x 
n 


+ x 

n 


+ hO)+mgwx -hkx 
n n n n n n n 


m n g0 = 0 


From these equations, the model force in the x-direction and the moment about 
the y-axis can be found (see Table 4-39). 
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Annular Tank | Spring-Mass Model 

Excitation: Harmonic Translation and/or Pitching 

Figure 4-3CK Equivalent Mechanical Model 
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Table 4-40. Model Elements 


Annular Tank 


Spring-Mass Model 


Excitation: Harmonic Translation or Pitching 


2 b 

Natural Frequency a) = ■“ £ tanh < 

M J nan n 


Spring Constant 


k = m to 
n n n 


(fig. 4-31) 


Damping Coefficient 
of Slosh Mass 


c = m g w 
n n °n n 


(fig. 4-32) 


Damping Coefficient 
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Ratio of 
Slosh Mass 
to Fluid Mass 


c 1 + 
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n __ n L n 1 n J n 


(1-k^) K 


(fig. 4-33) 
(fig. 4-34) 


Ratio of 
Fixed Mass 
to Fluid Mass 


Ratio of Slosh 
Mass Coordinate 
to Fluid Depth 


Ratio of Fixed 
Mass Coordinate 
to Fluid Depth 


Moment of Inertia 
of Solidified Fluid 


Moment of Inertia 

of Disc (c, = o) 
a 


Moment of Inertia 

of Disc (c. f o) 
a 


Moment of Inertia 
of Fixed Mass 
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,2 V 

ni L h 


(fig. 4-35) 


(fig. 4-36) 


(fig. 4 - 37 ) 
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DAMPING CONSTANTS 
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Ratio of Fluid Depth to Tank Radius 





Ratio of FI 
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4-35. Model Element Graph 
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Figure 4-36. Model Element Graph 




4-102 


Ratio of Fluid Depth to Tank Radius 
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Figure 4-38, Model Element Graph 
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Ratio of Fluid Depth to Tank Radius 



Table 4-41. Model Analysis 


Annular Tank 


Pendulum Model 


Excitation: Harmonic Translation and/or Pitching 


Figure 4-40shows a diagram of the pendulum model used in representing the 
dynamic response of a liquid in an annular tank when subjected to harmonic 
translation in the x-direction and/or pitching about the y-axis. 


Coordinate System: 

The origin is located at the center of gravity of the undisturbed liquid 


Model Description: 

The components of the system are as follows: 

1. A fixed mass M having a moment of inertia I is rigidly connected to 
the tank and is located on the z-axis at a distance H below the 
coordinate origin. 

2. A set of movable masses is distributed along the z-axis when the 
tank is at rest. These modal masses are pendulums having massless 
lever arms of length attached to the z-axis at distances H n above 
the origin. They are constrained by dashpots having viscous damping 
coefficients c n to remain approximately in the xy-plane and to move 
approximately parallel to the x-axis. Angular displacements of the 
pendulum with respect to the tank (z-axis) are denoted by 

3. A massless disc having a moment of inertia 1^ is loacted at the 
coordinate origin. Its motion is confined to rotation about the 
y-axis and is cushioned by a dashpot having a viscous damping 
coefficient c^. The angular displacement of the disc relative to 
the tank is defined by i|>. 


Equations of Motion: 

The equations, obtained through Lagrange T s equations, are as follows: 
1. Force Equation: 

® r 

F = - M(x + H8) - 7 m x + L X + (H - L )0 

x n L n n nn 

n=l 


2 . 


Moment Equation : 

M = - (I + MH 2 )6 - I J (6 + *) + g 
y d 

LX + (H - L 00 
n n ^ » 



m L X 
n n n 


00 


l 

n=l 


m (H 
n n 


L ) 


n 
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Table 4-41. Model Analysis (continued) 


Annular Tank 


Pendulum Model 


Excitation: Harmonic Translation and/or Pitching 


Equations of Motion (continued): 
3 . Disc Equation : 

1 (0 + ^) -f c ^ = 0 

d d 


4. Slosh -Mass Equation: 

r .. • 

m x + L A + (H - L )0 + m g w L A - ni gQ = 0 

n L n n n n J n n n n n n 


From these equations, the model force in the x-direction and the moment about 
the y-axis can be found (see Table 4-42). 
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Annular Tank 


Pendulum Model 


Excitation: Harm to ic Translation and/or Pitching 


Figure 4-40 . Equivalent Mechanical Model 
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Excitation: Harmonic Translation and/or Pitching 
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Table 4-43 . Model Elements 


Annular Tank 


Pendulum Model 


Excitation: Harmonic Translation and/or Pitching 


Natural Frequency 


Pendulum Length 


Damping Coefficient 
of Slosh Mass 


Damping Coefficient 


Ratio of 
Slosh Mass 
to Fluid Mass 


Ratio of 
Fixed Mass 
to Fluid Mass 


J 2 ~ £ tanh < 

nan n 


L 

” “g 



m A [2/ir£ - kC.(o )| tanh < 

n nL n 1 n J n 


(l-k Z ) k 


, , iu 

~ = i - I -=• 
”L n-o "l 


Ratio of Pendulum 
Mass Coordinate 
to Fluid Depth 


Ratio of Fixed 
Mass Coordinate 
to Fluid Depth 


Moment of Inertia 
of Solidified Fluid 


III - L I 1 


(fig. 4-42) 



Moment of Inertia 

of Disc (c, = o) 
d 


Moment of Inertia 

of Disc (c, f o ) 
a 


Moment of Inertia 
of Fixed Mass 


(1-K ) « 



2 i 2 

nub 


,2 u 2 

m^h m L^ 


T d M ,H> 

u 2 " in. n } 
m L h *L 


l - 

rA> “l 
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HEIGHT RATIO FOR CYLINDRICAL TANK 
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Figure 4-41. Model Element Graph 





SUSPENSION POINT LOCATION 
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Figure 4-42. Model Element Graph 




4.7 Annular-Sector Tank 
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z 



Container : The tank is a sector of an annular tank having an inner radius 

r=% and an outer radius r=a and is filled with a liquid to a depth h. 

Coordinate System : The origin is located at the center of gravity of the 

undisturbed fluid which would be contained in an annular tank generatec 
by revolving the annular-sector tank about the z-axis. The x-axis 
must lie in the sector wall. 

Reference : (4) 

Comments : The results given in this section are not valid for 5 = — , — . 
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Table 4-44. Boundary Conditions, Velocity Potential, and Natural Frequency 
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Boundary Conditions, Velocity Potential, and Natural Fluency (continued) 



ij--i±5 


Table 4-45. Liquid Force and Moment Resultants 
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Table 4-45. Liquid Force and Moment Resultants (continued) 
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Table 4-46. Boundary Conditions, Velocity Potential, and Natural Frequency 
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Boundary Conditions, Velocity Potential, and Natural Frequency (continued) 
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where F are roots of A (f ) = J (f ) Y (k£ ) - J (kr ) Y (r 

mn m/2a 4 mn m/2a ^mn m/2a ^mn m/2a ^nm m/2a Sn 



Annular -Sector Tank | Excitation: Harmonic Pitchin 

iu)t iu)t 

Pitching about the y-axis, 0 = 0 e and pitching about the x-axis, x = X e 
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Table 4-47* Liquid Force and Moment Resultants (continued) 
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Table 4-48. Boundary Conditions, Velocity Potential, and Natural Frequency 
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Table 4-49. Liquid Force and Moment Resultants 
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Table 4-49. Liquid Force and Moment Resultants (continued) 
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Table 4-49* Liquid Force and Moment Resultants (continued) 
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4.8 Rectangular Tank 
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a 


Co ntainer : The tank is rectangular having base dimensions a and b in the 

x and y directions, respect ively, and is filled with a liquid to a depth 

h . 


Coordinate System ; The origin is located at the center of gravity of the 
undisturbed fluid. 

References : (76) 
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Table 4-50. Boundary Conditions, Velocity Potential, and Natural Frequency 
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(2n+l)(7r/a) tanh ( 2 n+l) — 



Liquid Force and Moment Resultants 
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Table 4-52. Model Analysis 


Rectangular Tank 


Spring-Mass Model 


Excitation: Harmonic Translation and/or Pitching 


Figure 4-43 shows a diagram of the spring-mass model used in represent- 
ing the dynamic response of a liquid in a rectangular tank when subjected to 
harmonic translation in the x-direction and/or pitching about the y-axis. 


Coordinate System: 

The origin is located at the center of gravity of the undisturbed 
liquid. 


Model Description : 

The components of the system are as follows: 

1. A fixed mass M having a moment of inertia I is rigidly connected 
to the tank and is located on the z-axis at a distance H below the 
coordinate origin. 

2. A set of movable masses m is distributed along the z-axis when 
the container is at rest at distances h above the origin. These 
modal masses are constrained by springs^aving spring stiffness 
coefficients k to remain in the xy-plane and to move only in a 
direction parallel to the x-axis. 


Equations of Motion: 

The equations, obtained from either an equilibrium or energy formulation, 
are as fo 1 lows : 

1. Slosh Mass Equation for Translation: 

•• 2 
mx+kx^mxo) sinwt 
n n no 

2. Slosh Mass Equation for Pitching: 


mx+kx = (mg+hk')0 sinwt 
n n n n n o 

From these equations, the model force in the x-direction and the moment about 
the y-axis can be found (see Table 4-53). 


* 
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Rectangular Tank | Spring -Mass Model 

Excitation: Harmonic Translation and/or Pitching 

Figure 4-43. Equivalent Mechanical Model 


4-132 






Excitation: Harmonic Translation and/or Pitching 
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Table 4-54. Model Elements 


Rectangular Tank 


Spring-Mass Model 


Excitation: Harmonic Translation and/or Pitchin- 


Natural Frequency 


2 = S. 


(2n+l)*rr tanh £(2irH)iTr^ 


Spring Constant 


Ratio or Spring 
Constant to Fluid 
Weight 


Ratio of Slosh 
Mass to Fluid Mass 


k = m w 
n n n 


hk^ 8 tanh^ £(2n+l )irr ^ j 


m L 8 [(2n+l)7t] 2 


8 tanh ^(2n+l Vr ^ j 
m L £(2n+l )i T ] 3 r 1 


(fig. 4-44) 


(fig. 4-46) 


— = i - z — 

m u -m 


Ratio of Fixed 
Mass to Fluid Mass 


Ratio of Slosh Mass 
Coordinate to 
Fluid Depth 


Ratio of Fixed 
Mass Coordinate 
to Fluid Depth 


Moment of Inertia 
of Solidified Fluid 


2 2 

Effective Moment of m L^ 

Inertia of Fluid 


m T 

n^ L 


h n | ^ tanh 2n+l )rr r ^ / 2 j 
h 2 (2n+l)nr 1 /2 


IhI = ^ y 

h M *• 



(fig. 4-46) 


(fig. 4-45) 


(fig. 4-45) 



_ E . + Z68_ 

2 2, . . 2, .. 5 

L r + 1 (r^+ l)* Tj 

tanh[(2n+l > 1 ^/ 2 ] ' 

n=o ( 2n+l )"* 


(fig. 4-47) 


Moment of Inertia 
of Fixed Mass 


2 = ~ 2 “ m~ ( l )2 - l ( ^ )2 (fi8 ‘ 4 " 47) 

l i l m T h L m T h 


m^h m^h 


n=o L 
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NON- 


k 0 = SPRING CONSTANT ASSOCIATED WITH FUNDAMENTAL 

k , = SPRING CONSTANT ASSOCIATED WITH FIRST ODD 
HARMONIC 

k 2 = SPRi NG CONSTANT ASSOCIATED WITH SECOND 

ODD HARMONIC 

W L = WEIGHT OF FUEL I 



BIE3 


0.5 1.0 1.5 2.0 2.5 

TANK ASPECT RATIO r, 

Non-Dimensional Spring Constants 
vs Tank Aspect Ratio. 


Figure 4-44. Model Element Graph 
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RATIO OF ARMS TO DEPTH 



H * VERTICAL LOCATION OF FIXED MASS 

h 0 = VERTICAL LOCATION OF FUNDAMENTAL 
MOVING MASS 

VERTICAL LOCATION OF FIRST ODD 
HARMONIC MOVING MASS 




0.5 i.O 1.5 2.0 2.5 

TANK ASPECT RATIO r, 

Ratios of the Vertical Locations or Arms 
of the Fixed Mass and First Two Slosh 
Masses to Fluid Depth vs Tank Aspect Ratio. 


Figure 4-45. Model Element Graph 
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h~i 




Ratio of Fixed Man end First TWo Slosh 
Masses to Fluid Mass vs Tank Aspect Ratio. 


Figure 4-46. Model Element Graph 
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Figure 4-47. Model Element Graph 



Table 4-55. Boundary Conditions, Velocity Potential, and Natural Frequency 
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Table 4-56. Liquid Moment Resultant 



k-lko 



Table 4-57. Model Analysis 


Rectangular Tank 


Excitation : 


Torsional Spring-Mass Model 


Harmonic Roll 


Figure 4-48 shows a diagram of the torsional spring-mass model used in 
representing the dynamic response of a liquid in a rectangular tank when sub- 
jected to harmonic roll about the z-axis. 


Coordinate System: 

The origin is located at the center of gravity of the undisturbed liquid. 


Model Description : 

The components of the system are as follows: 

1. A fixed mass M having a moment of inertia I is rigidly connected to 
the tank and is located on the z-axis at a distance H below the 
coordinate origin. 

2. A set of movable masses m is constrained by torsional springs 

having stiffness coefficients k to remain in planes parallel to the 

xy-plane and at distances h a^ove it. 

J * mn 


Equation of Motion: 

The equation, obtained from either an equilibrium or energy formulation, 
is as follows : 

«« 

I _♦ + = k sinwt 

mn mn mn o 

From this equation, the moment about the z-axis can be found (see Table 4-58). 
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Figure 4-48. Equivalent Mechanical Model 
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Table 4-58. Model Moment Resultant 



4 - 14-3 


Table 4-59. Model Elements 


Rectangular Tank Torsional Spring-Mass Model 

Excitation: Harmonic Roll 

Natural Frequency 

0 ) 2 Vr 2 (2m+1) 2 + (2n+l) 2 

mn h 3 

tanh [ifr j Vr 2 (2m+1) 2 + (2n+l) 2 J 

Torsional Spring 
Constant 

2 

k = x co 

mn mn mn 

Moment of Inertia 
of Solidified Fluid 

_ m 2 , , 2 

X s " U a + b 

Effective Moment of 
Inertia of Fluid 

i 4 768 y tanh[(2n+l)irr 3 /2] 

T 1 ” o 5 o Z 5 (f ig*4- 

i s r 3 +l IT r 3 (r 3 +l) n=^o (2n+l) 49,4-50) 

Moment of Inertia 
of Slosh Mass 

I 768 [r?(2m+l ) 2 - (2n+l) 2 ! 2 

mn L J J , 

T S TT 6 (HT 2 )(2m+l) 4 (2n+l) 4 [r 2 (2m+l) 2 + (2n+l) 2 ] 

tanh[u ri Vr 2 (2m+l) 2 + (2n+l) 2 ] ^ 

Vr 2 (2m+1) 2 + (2n+l) 2 4 " 50) 

Moment of Inertia 
of Fixed Mass 

_ 00 CO J 

I _ I y y mn 

I I " 1 L I 

! s s m— o n—o s 
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MOMENT OF INERTIA RATIOS 



TANK ASPECT RATIO r. 


Roll Moment of Inertia Ratios for 
Shallow Tank vs Tank Aspect Ratio. 


Figure 4-49. Model Element Graph 
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Figure 4-50. Model Element Graph 
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5.1.1 Introduction . Sloshing becomes detrimental when the wave 


tmplitudes of tanked propellants, agitated by oscillations, approach 
breaking height. This sloshing of the liquid exerts considerable forces 
on the missile structure, affects the controls and renders any liquid 
level measuring device inaccurate. This means that a sloshing damping 
device has to be provided. The term damping usually refers to the 
checking of a motion due to resistance, as by friction or any other 
similar cause. It is of especial significance in connection with the 
diminishing amplitude of an oscillation, as that of a liquid swinging 
Inside a moving container. Unless energy is supplied during each cycle, 
the amplitude of the fluid motion falls off at each successive oscillation 
by an amount commonly expressed in terms of the decrement, or damping 
factor, which is the ratio of any one amplitude to that next succeeding 
it in the same sense and direction. 

Many of the force and moment resultants given in Chapter IV require 
a damping factor before they can be evaluated. The accuracy of these 
equations and consequently their usefulness is therefore very much dependent 
on the value chosen for the damping parameter. When a large propellant 
tank has no damping devices the resistance to fluid motion at the container 
walls is generally accepted to be nil. That is, the fluid damping is 
negligible. However, when damping systems are introduced, the resulting 
damping factor depends on many variables and is thereby extremely difficult 
to obtain. For some damping devices, such as fixed ring baffles, a sound 
analytical approach can be effectively used while in most cases an experi- 
mental approach is the only way of obtaining a meaningful damping factor. 
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5.1.2 Design Requirements of an Anti-Slosh Device . Studies 
(40), (42) show that the depth below the undisturbed free surface at 
which the liquid no longer participates in surface oscillations is 
approximately equal to one-fourth the tank diameter. Therefore a damping 
device need only extend to such a depth. The liquid level, however, 
changes as the container empties thereby requiring either a damping device 
capable of moving down at a rate equal to the draining rate or a damping 
device which is fixed to the wall but extends to the bottom of the con- 
tainer. The first stipulation suggests some type of floating device, 
whereas the second suggests either longitudinal partitions or an evenly 
spaced baffle system covering the inner periphery of the tank and distri- 
buted throughout the depth of the liquid. 

The amount of damping needed to stabilize a missile is deter- 
mined from analog or digital computer studies of the over-all dynamic 
characteristics of the vehicle. In some cases, stability studies might 
indicate no need for anti-slosh baffles in a particular tank. Generally, 
however, the designer is faced with the problem of optimizing a damping 
system by integrating into his design the requirements that follow: 

11 Produce a high damping effect. 

2. Pay a minimum weight penalty. 

3# Utilize minimum space. 

4. Absorb the liquid forces and moments or transfer them uniformly 
to the tank structure, thereby avoiding points of high stress 
concentration. 


5-3 



5* Function throughout environmental changes, such as temperature. 

6. Not interfere with other operations, such as the emptying of the 
container or the measurement of the liquid level. 

7. Be easy to assemble. 

8. Not cause damage to the tank or other built-in equipment during 
transportation of the missile. 

9. Not interfere with cleaning operations. 
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Ring Baffles 



where C = amplitude of liquid measured at the tank wall 
h = depth of liquid 

d = depth of ring below undisturbed free surface 
a = radius of right circular cylinder 
w = ring width 

References : (32, 35, and 41)* 

Ring baffles are most effective against lateral slosh which has large 
vertical components of fluid velocity near the tank wall. 

Investigations by Silveira, Stephens and Leonard (35) indicate that-: 
1. For a given baffle depth: 

(a) The damping factor increases as the baffle width increases 
(fig* 5-10). 

(b) The frequency increases as the baffle width decreases 
(fig. 5-15). 
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Ring Baffles (continued) 


2. Based on the total surface of a baffle, the highest mean damping 
factor appears to be afforded by ring baffles as compared to the 
baffle types depicted in fig. (5-77). 

Investigations by Garza and Abramson (32) indicate that: 

1. For a given baffle depth: 

(a) With the hole size constant: 

(1) The frequency increases for increasing percent per- 
foration (fig. 5-2) . 

(2) The damping factor increases for decreasing percent 
perforation (fig. 5-8). 

(b) With the percent perforation constant: 

(1) The frequency increases for increasing hole size 
(fig. 5-3). 

(2) The damping factor increases for decreasing hole size 
(fig. 5-7). 

2. Damping produced by perforated baffles is consistently lower 
than that produced by a solid baffle (fig. 5-7 and 5-8). 
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5.1;3, Ring Baffles * Because of the complexity of the boundary con- 
ditions, a theoretical approach for predicting the damping of a device is 
extremely limited. In particular, a fixed ring baffle is the only damping 
device for which an analytical approach exists. For this reason, it is 
necessary to resort to an experimental approach consisting of ground 
oscillation testing of a full scale or model tark containing the damping 
device to be evaluated. If the size of the model is small as compared to 
the prototype, the use of similitude theory is advisable in order to obtain 
reasonable similarity of fluid motion. Since it is generally not possible 
to model equal values of all dimensionless parameters affecting the fluid 
motion, an appropriate nondimensiona 1 modeling parameter must be chosen. 
Viscous forces and surface tension are small and can therefore be neglected 
whereas, forces resulting from both the velocity of the fluid and the 
inertial effects of acceleration cannot. Acceleration is of especial 
importance because of the three following effects on slosh phenomena: 

1. During ground tests, the eigen frequencies remain constant, 
whereas, during flight conditions they vary as the square 
root of the longitudinal acceleration. 

2. The wave amplitudes are approximately inversely proportional 
to the longitudinal acceleration. 

3. The damping factor is proportional to the square root of the 
liquid amplitude. 

4. The force and moment resultants are proportional to g. 

Therefore, Froude T s number, being the ratio of fluid velocity forces 
to inertia forces, is selected for the modeling parameter. 
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5.1.3. 1 Ring Baffles: Experimental . Each of Articles 5. 1.3. 1.1 


through 5. 1.3. 1.3 presents an analysis concerned with a particular reference 
from which experimental data in the form of graphs is presented. The analyses 
are composed of two parts; the first being concerned with experimental 
procedure and the second with interpretation of results. This information 
is meant to aid the designer in both evaluating data of a particular source 
and comparing data of different sources. Also, this information should be 
of use in similitude applications since most of the data is obtained from 
model studies. 


5. 1.3. 1.1 Investigations by Garza and Abramson (32) are 
outlined below. 

Fixed Ring Baffles 

1. Baffle thickness: 0.018" to 0.030". 

2. Ratio of baffle width to tank radius, ^ = 0.157. 

R 

3. Perforation: 

(a) The hole size was 0.079" while the percent perforation 
was varied. 

(b) The percent perforation was 30% while the hole size 
was varied. 

4. Attachment to container: The baffles were secured to the 

tank by four 1/8" x 1/3" steel strips attached at 90° apart 
around the outer edge of the ring. They were supported 
through angle iron brackets bolted to the upper flange of 
the tank. The outer edge of the ring was turned down, 
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giving the ring an L type cross section which stiffened the 
ring considerably, but had little effect on damping (see 
Figure 5-1). 

Container, Excitation and Instrumentation 

The container was a rigid-wall circular cylindrical tank sup- 
ported by four dynamometers . All tests were conducted for three amplitudes 
of translational excitation. Because of the rather significant effects of 

excitation amplitude, measurements were made for various values ranging from 
x 

0.00184 <_ jc 0.00823, and then all data was presented in terms of RMS 
values. Table 5-1 gives the resonant frequencies for each of the various 
tests conducted. 

Damping Factor 

Experimental values of the damping factor were obtained from 
resonant peaks of experimental force response curves. Theoretical values 
of Y g were calculated from Miles f equation (Article 5. 1.3.2. 1) using liquid 
surface amplitudes measured at the tank walls. Table 5-2 gives the 

damping values for each of the various tests conducted. 

5. 1.3.2 Investigations by Silveira, Stephens, and Leonard (35) 
are outlined below. 

Baffles (Figure 5-76) 

1. Baffle thickness 

(a) 12 inch tank: conic section, 1 / 16 ” Plexiglas; all others, 

1/8" Plexiglas. 
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(b) 30 inch tank: fixed ring, 1/4" Plexiglas. 

2. Ratio of baffle width to tank radius 

(a) 12 inch tank: see Figure 5-77. 

(b) 30 inch tank: W/R = 0.076. 

3. Perforation 

(a) Hole size, 1/8 n diameter. 

(b) Percent perforation, 50%. 

Container Excitation and Instrumentation 

Two containers were used: a 12" diameter tank of 1/8" Plexiglas 

and a 30" aluminum tank of 0.016" walls and 1/2" base. A paddle was used 

to excite the liquid in the fundamental mode and then removed when the 

amplitude was sufficient. Liquid response in the 12" tank was sensed by 
strain gages mounted on the torsion bars which are located between the 
base platform and the support ring, whereas the response in the 30" tank 
was measured by a load cell which replaced one of the three platform supports. 
The output signals from the strain gauges and the load cell were amplified 
and then fed into a dampometer to measure the damping and frequency of 
the liquid motion. 

Damping Factor 

The rate of decay of the moment resulting from the damped liquid 
oscillations was measured and the damping factor <5, representing the decay 
of oscillation, was defined as 
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6 


M 

„ lo8 e M n 


where n is the number of cycles over which the decay was measured, is 
the magnitude of the initial moment and is the magnitude of the moment 
after n cycles. 

General Comments 

1. No graph is given for cruciform baffles since, for all but 

shallow depths, damping is independent of the fluid depth. 
Test results for the 90° and 45° positions (see Figure 5-77) 
are as follows: for ^ = 0.169, 6^ = 0.072 and 6^ = 0.070 

and for j- = 0.337, 6^ = 0.156 and 6^ = 0.142. 

2. All data points presented represent the average of five or 
more measured values for the given condition. 

3. For high values of baffle depth, the curves in Figures 5-9, 
5-10, 5-11, 5-66, 5-67, and 5-69 approach the damping values 
for which no baffle is present* 

4. In Figure 5-77, the mean damping factor for a particular 

baffle was obtained in the following manner: from the graph 

the damping 

relating/factor to baffle depth, the area under the curve 
d d 

between 0.084 above and 0.084 — below the maximum damp- 
K , K 

a 

ing was divided by 0.168 . 

K 

5. The curves in Figures 5-12 and 5-14 are given by Miles' 
equation as presented in reference 33. It is 
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6 = 


- -4.6- 

5.66ire R 


R 2 - (R-W) 2 


>2 


2 (r)^ 


If W is neglected, then the damping factor becomes: 


A . -4.6^ n 1 w\2 

o = 5.66*e R <* \R~' 


which is identical to Miles' equation except that the 
constant multiplier is 5.66 if instead of 2.83. 

6. For high values of baffle depth, the curves in Figures 5-15, 
5-16, 5-64, 5-65, and 5-68 approach, asymptotically, the first 
resonant frequency for a circular cylindrical tank con- 
taining no baffles. 
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Figure 5.1 Details of battle support arrangements ana tank configuration. (32) 
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TABLE 5.1 


Liquid Resonant Frequency. (32) 
W/R « 0. 157 


% Open Area 

0 

8% 

16% 

23% 

30% 

23% 

30% 

30% 

Hole Diameter 

(in. ) 

0 

0. 079 

0. 079 

0. 079 

0. 079 

0. 040 

0. 040 

0. 02 

Ring Depth 
d s /R 



(jJ Z d/g 

(Xo/d* 

0.00184) 




0 

4.68 

4.40 

4. 25 

3. 57 

4. 1 

4. 11 

4. 14 

4. 33 

. 025 

4. 19 

3. 19 

3.47 

3. 30 

3. 88 

3. 20 

3.43 

3. 26 

. 050 

2. 92 

3. 17 

3. 31 

3. 34 

3. 52 

3. 36 

3. 34 

3.22 

. 075 

3. 01 

2.90 

3.27 

3. 32 

3.47 

3.26 

3. 34 

3. 32 

. 100 

3. 06 

2. 96 

3. 31 

3. 35 

3. 51 

3. 32 

3.41 

3. 34 

. 125 

2.97 

3. 19 

3.40 

3. 36 

3.49 

3. 31 

3.43 

3. 35 

. 175 

3. 18 

3. 24 

3. 35 

3. 32 

3. 54 

3. 38 

3.47 

3.41 

. 250 

3.40 

3. 36 

3.40 

3.45 

3. 56 

3. 52 

3. 54 

3.42 

. 375 

3.48 

3.49 

3. 54 

3. 58 

3. 57 

3. 52 

3. 56 

3. 54 

.450 

3.47 

3. 54 

3.62 

3. 58 

3. 58 

3. 59 

3. 62 

3. 54 




C0 2 d/g 

(X 0 /d» 

0.00417) 




0 

4.43 

3.84 

3.76 

3.73 

3.65 

3. 81 

4. 16 

3 87 

. 025 

4. 27 

3. 64 

3.79 

3.62 

3. 62 

3.60 

3. 54 

3. 90 

. 050 

3.71 

3.29 

3.40 

3. 54 

3. 50 

3.44 

3.66 

3. 55 

. 075 

2.99 

3. 14 

3.28 

3. 32 

3.46 

3. 26 

3.45 

3. 31 

. 100 

2. 94 

3. 15 

3.28 

3. 35 

3.45 

3.28 

3.38 

3. 30 

. 125 

3. 09 

3.20 

3. 22 

3. 35 

3.45 

3.26 

3.40 

3. 31 

. 175 

3. 09 

3. 24 

3. 26 

3. 38 

3. 51 

3.26 

3.45 

3. 34 

.250 

3. 18 

3. 34 

3. 38 

3.43 

3.49 

3. 39 

3.45 

3.40 

. 375 

3. 33 

3.46 

3.40 

3.43 

3.42 

3. 52 

3.45 

3. 44 




U) 2 d/g 

(Xo/d - 

0.00833) 




0 

3. 96 

3. 84 

3.73 

3. 54 

3. 52- 

3. 58 

3. 28 

3.49 

. 025 

3.77 

3.77 

3.75 

3.49 

3.50 

3. 54 

3.41 

3. 54 

. 050 

3. 58 

3. 58 

3. 51 

3.49 

3.47 

3.47 

3.40 

3.49 

. 075 

3. 18 

3. 24 

3. 28 

3.40 

3.40 

3.40 

3.41 

3.40 

. 100 

2.97 

3.22 

3. 28 

3. 31 

3. 38 

3.24 

3. 34 

3. 36 

. 125 

3. 00 

3. 06 

3. 32 

3.24 

3. 31 

3.28 

3. 30 

3. 28 

. 175 

3. 07 

3. 13 

3.28 

3.31 

3. 31 

3. 28 

3. 30 

3. 32 

.250 

3. 16 

3.28 

3.28 

3. 31 

3.28 

3.22 

3. 32 

3. 32 

. 375 

3. 30 

3. 32 

3.31 


3.24 

3. 30 

3. 36 

3. 32 
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TABLE 


5.2 


Ring Damping. (32) 
W/R * 0. 157 


% Open Area 


00 

16% 

23% 

30% 

23% 

30% 

30% 

Hole Diameter 

fin. 1 

0 

0. 079 

0. 079 

0. 079 

0. 079 

0. 040 

0. 040 

0. 020 

Ring Depth 
ds/R 


Damping Ratio, 

X 0 /d« 

0. 00184 




0 

. 060 

. 087 

. 084 

. 108 

. 087 


. no 

. 103 

. 025 

.210 

. 132 

. 084 

. 106 

. 087 

. 114 

. 112 

. 116 

. 050 

. 137 

. 147 

. 071 

. 073 

. 055 

. 086 

. 086 

. 103 

. 075 

.090 

. 1145 

. 070 

.065 

. 050 

. 072 

. 082 

. 084 

. 100 

.070 

. 0885 

. 069 

.059 

. 045 

. 046 

. 050 

. 057 

. 125 

. 073 

.067 

. 056 

. 049 

. 045 

. 059 

. 045 

. 053 

. 175 

. 060 

. 054 

. 056 

. 045 

. 040 

. 054 

. 033 

. 049 

.250 

.035 

. 040 

. 035 

. 047 

. 031 

. 032 

. 032 

.042 

. 375 

. 027 

. 025 

. 026 

. 021 

. 022 

. 023 

. 022 

. 029 

.450 

. 025 

. 022 



. 017 

. 023 

.015 

. 023 



Damping Ration 

X Q /d ■ 

0.00417 




0 

. 101 

. 127 

. 101 

.098 

. 082 

. 094 

. 102 

. 086 

. 025 

. 121 

. 126 

. 122 

. 108 

. 093 

. no 

. 114 

. 109 

. 050 

. 173 

. 134 

. 119 

. 112 

. 089 

. 109 

. 106 

. 103 

.075 

. 160 

. 112 

. 108 

. 084 

. 077 

. 104 

. 077 

. 089 

. 100 

. 110 

.091 

. 090 

. 075 

. 064 

. 085 

. 068 

. 085 

. 125 

. 077 

. 102 

. 077 

. 064 

. 047 

. 077 

. 064 

. 078 

. 175 

. 072 

. 070 

. 072 

. 061 

. 050 

• 066 

. 054 

. 070 

. 250 

. 051 

. 061 

. 051 

. 049 

. 047 

. 051 

. 049 

. 055 

. 375 

. 039 

. 044 

. 038 

. 035 

. 040 

. 039 

. 032 

. 037 



Damping Ratio, 

X 0 /d- 

0.00833 




0 

. 117 

. 122 

.121 

. Ill 

• .095 

. no 

. 091 

. 114 

. 025 

. 146 

. 126 

. 124 

. 117 

. 087 

. 122 

. 106 

. 112 

. 050 

. 177 

. 134 

. 124 

. 109 

. 098 

. 126 

. 098 

. 118 

. 075 

. 158 

. 123 

. 114 

.099 

. 084 

. 104 

. 086 

. 120 

. 100 

. 142 

. 115 

. 102 

. 093 

. 082 

. 093 

. 072 

.093 

. 125 

. 116 

. 106 

.093 

. 088 

. 078 

. 093 

. 082 

. . 088 

. 175 

. 096 

. 096 

. 081 

. 078 

. 072 

. 079 

. 074 

. 072 

.250 

. 082 

.066 

.067 

. 061 

. 057 

. 067 

. 058 

. 058 

. 375 

. 066 

.048 

.046 


. 026 

. 046 

. 038 

. 042 
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RING 
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RMS ( Xo/d « 0.00184 to 0.00833) 
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BAFFLE DEPTH , d # /R 


Figure 5*2 Effect of percent perforation on liquid resonant 
frequencies as a function of baffle depth. (32) 











LIQUID RESONANT FREQUENCY PARAMETER, 
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- 
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Figure 5.3 Effect of perforation hole size on liquid resonant 
frequencies as a function of baffle depth. (32) 





BAFFLE DEPTH, d,/R 
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DAMPING RATIO 



Figure 5-5 Comparison of theory and experiment for damping provided by a 
flat solid ring baffle as a function of b a ffl e depth. (32) 
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DAMPING RATIO, 


0.32 


I SYMBOL 

0.28 1 G 


0.24 h 


Xo/d 

0.00184 

0.00417 

0.00833 

RMS (Xo/d - 
0.00184 to 
0.00833) 



h*2R 


0.20 


0.16 


W/R* 0.157 





0.08 


0.04 






0.08 0.16 0.24 0.32 

BAFFLE DEPTH, d $ /R 


0.40 


Figure 5.6 Effect of excitation amplitude on damping effectiveness of a 
solid flat ring baffle as a function of baffle depth. (32) 


5-20 


DAMPING RATIO 
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0.24 
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Figure 5.7 Effect of perforation hole size on damping 

effectiveness as a function of baffle depth. (32) 
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BAFFLE DEPTH, d,/R 


Figure 5.8 Effect or percent perforation on damping 

effectiveness as a function of baffle depth. (32) 
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Baffle location, -^s 


Variation of damping factor with baffle location for fixed 

ring baffle. R = 6 inches; - = 2. (35) 

R 


Figure 5-10 













Variation of damping factor with baffle location in the 
12- and 30-inch-diameter tanks. (35) 


Figure 5-12 
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Figure 5-13 







Baffle location, -^6 


Variation of frequency with baffle location for fixed-ring 

baffle. (35) 

Figure 5-15 
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Frequency, cps 



Variation of frequency with baffle location for ring-with- 
radial-clearance baffle. (35) 


Figure 5-16 
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5. 1.3.2 Ring Baffles: Analytical . Damping requirements are 

determined from stability studies of the vehicle in question. The designer 
must then determine the baffle width and spacing needed to yield a damping 
factor in the required range. 

5. 1.3. 2.1 Miles 1 Analysis (28) . According to Miles, the damp- 
ing ratio of an annular ring baffle is 

4.6d 

y = 2.83 a 3/2 e 3 (^/a)^ 2 , a«l and <^<<3 (5-1) 


where 

2 

a7r a = ring baffle area 

a = ratio of ring baffle area to cross-sectional area of tank 
d = depth of ring baffle below undisturbed free liquid surface 
a = radius of circular cylindrical tank 

^ = amplitude of liquid oscillations measured at the tank wall 
from undisturbed free liquid surface. 

Within the restrictions on a and this damping ratio relation was exper- 
imentally verified ; however, Miles suggested that additional experimental 
confirmation would be desirable, both in order to establish the limits on 
a and C^/a and to determine a more accurate value of the constant multiplier. 
Results are shown in Figures 5-17 to 5-29. A nomograph is shown on page 5-32. 
A line drawn between values of C^/a and d/a determines a pivot point on the 
K-axis. A line drawn from this pivot point through a value of 7 gives the 
corresponding value of W/a ? 
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MILES' DAMPING RATIO 
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DAMPING FACTOR VERSUS 
LIQUID AMPLITUDE £ w “0 .< 



Figure 


5-3 



% 



Figure 5-19 
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DAMPING FACTOR VERSUS BAFFLE DEPTH FOR WIDTH 
W* 0.025 WITH LIQUID AMPLITUDE AS A PARAMETER 
(MILES FORMULA) 



Figure 5-22 
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OAMPMQ FACTOR VERSU8 BAFFLE DEPTH FOR WDTH 
W-0.050 WITH LIQUID AMPLITUDE AS A PARAMETER 
(MILES FORMULA) 



Figure 5-23 
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DAMPING FACTOR VERSUS BAFFLE DEPTH FOR WIDTH 
WO.IOQ WITH LIQUID AMPLITUDE AS A PARAMETER 
(MILES FORMULA) 



Fieure 5-24 
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Figure 5-25 

5-4-1 


DAMPING FACTOR VERSUS SURFACE AMPLITUDE FOR BAFFLE 

WIDTH W« ,05a WITH BAFFLE LOCATION AS A PARAMETER 

(MILES FORMULA) 



5-42 






DAMPING FACTOR VERSUS SURFACE AMPLITUDE FOR BAFFLE 

WIDTH W-.l5a WITH BAFFLE LOCATION AS A PARAMETER 
(MILES FORMULA) 



Figure 5-27 



DAMPING FACTOR VERSUS BAFFLE WIDTH FOR 
£ w *.IOo WITH LOCATION OF BAFFLE AS A 

PARAMETER (MILES FORMULA) 
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Figure 5-28 
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DAMPING FACTOR VERSUS BAFFLE WIDTH FOR 

£ * ,30a WITH LOCATION OF BAFFLE AS A 
*W 



Figure 5-29 
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T>I <0 


5. 1.3. 2. 2 O'Neill's Extension (84). Within the scatter 


of experimental data, O'Neill showed that either of 

_ p 

can be used for reasonable values of — and F = 

a Pga 

£ o and for ring widths corresponding to a < 0.25 


the following equations 
y, for any ring submergence 


V = 2.83 e - 4 ’ 60d / a 

« 3 ' 2 

y = 2.16 e - 4 - 6(M /a 

a 3 ' 2 



(Miles) 


(5-2) 


Equation (5- 2) gives the damping factor in terms of the lateral liquid force F 
which is easier to measure than the wave amplitude Cj* 
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5*1*3. 2. 3 Bauer T s Extension (31) . During sloshing, the 
baffles and the oscillating fluid are not in continuous contact with each 
other* There are parts of the oscillation cycle during which the fluid does 
not uniformly wet the entire surface of the baffles. For this reason, Miles 
(28) expression, Equation ( 5-1), gives a higher damping effect than that 
experimentally observed. To account for the deviation, an effective baffle 
area concept was herewith successfully employed. This effective baffle area 
is a function of the width w of the baffle, its location d below the liquid 
free surface, and the maximum amplitude of sloshing The theoretical 

formula obtained agrees with experimental results. 

In Equation ( 5-1 ), the baffle area ratio blocking the 
cross sectional area is defined as 



if the baffle is completely submerged during a slosh cycle. If the baffle 
is completely out of the liquid during a certain time of the slosh cycle, 
the effective baffle area is 


d/a x , 2 

V a * 


C 1 /a + 
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5. 1.3. 2. 3 


Bauer f s Extension (31) (continued) 


^ + 7? 

d/a 

1 - w/a 

d 

i / (c./a) 2 - - 

T (C^/a) (1 - w/a) 


2 

7T 


f (1 - f) 

a a 

Cj/a 


In 


- ± a - ?) : 


arc sin 


d/a 


(1 - f); x /a 


+ \ (1 


Cj/a 


(5-3) 


d **1 

If the baffle is completely submerged, i.e., — > — , then the value 

a — a 

of a is equal to the underlined term (Equation 5-3). For the case of the 
baffle being out of the liquid during the slosh cycle in such a fashion 
that the inner rim of the ring baffle is not out of the fluid at any time, 
i.e. , 



a 



— a 


< 



a 


the underlined dotted terms are added to obtain the effective baffle ar i. 
For the case in which part of the baffle gets out of the liquid, i.e.. 


d 

a 



the total formula applies. Thus, Miles Equation for a single flat ring 


baffle becomes 


5. 1.3. 2.3 Bauer’s Extension (31). (continued) 


4.6d , 

Y = 2.83 (a) 3/ e" a (q/a) 2 (5-4) 

For a system of ring baffles, the damping factor can be obtained by super- 
imposing the contribution of each baffle. The nth baffle at a location 
d + (n-l)D below the free fluid surface exhibits an effective baffle area of 



tt q/a 
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5. 1.3. 2. 3 Bauer’s Extension (31) (continued) 
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5. 1.3. 2. 3 Bauer 1 s Extension (31) (continued) 


4 +(„-!) £<-!•(! . 2 ) 

a a a a 


the complete formula applies. 


For a baffle located above the undisturbed liquid surface, 
the same formula can be applied approximately with a slight modification. 

If the baffle is completely out of the liquid, i.e*, 



a 


< 


* 

d_ 

a 




where d = (D - d) is the distance of the baffle above the undisturbed 
liquid surface, the baffle area subjected to the fluid is a =0. The 
value D is the distance between baffles. For 



a a i w 
a 


and 



a 


* 

> d - J a . 



a 


9 


i.e., only a part of the baffle is subjected to the liquid during a slosh 
cycle, then the effective baffle area contributing to the damping is 
a = a - a. For the mth baffle, this is 
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5. 1.3. 2.3 Bauer f s Extension (31) (continued) 


If the value 


m 


D 

a 



9 


then = 0, while for the case where a partial part of the baffle is 
only subjected to the liquid during one slosh cycle, the effective area 
contributing to damping is presented by the underlined part of the 
formula. This means that this part is used for 


fi a .”)<£. i < ii 

a a — a a — a 


D d l/i w \ 
m — - — < — (1 - 
a a — a a 


the total formula is applied for the effective area. Thus the total damping 
of n baffles of width w submerged into the undisturbed liquid and m baffles 
of the same width outside of the undisturbed fluid, all of which are apart 
from each other by the value D, is given by 



n=l m=l 


(5-8) 


The result of this investigation can be seen in Figures 5-31 to 5-63. 
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Figure 5-30 

5-54 




U- i 




Figure 5-31 
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ANNULAR RING BAFFLE PARTI A l LY OUT OF LIQUID 
DURING SLOSH CYCLE 








EFFECTIVE BAFFLE AREA VERSUS BAFFLE DEPTH 

d FOR LIQUID AMPLITUDE C £ 0.5o AND VARIOUS 
BAFFLE WIDTH W 
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rigure 5-33 
















Figure 5-36 
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DAMPING FACTOR VERSUS BAFFLE DEPTH FOR 
LIQUID AMPLITUDE ^0 050 WITH THE BAFFLE 
WIDTH AS A PARAMETER (MILES FORMULA) 



Figure 5-37 
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DAMPING FACTOR VERSUS BAFFLE DEPTH FOR 
LIQUID AMPLITUDE C «O.IOa WITH THE BAFFLE 

WIDTH AS A PARAMETER (MILES FORMULA) 
REVISED 



Figure 5-38 
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DAMPING FACTOR VERSUS BAFFLE DEPTH FOR 
LIQUID AMPLITUDE 0.20a WITH THE BAFFLE 

w:dth as a parameter (miles formula) 



Figure 5- 39 
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DAMPING FACTOR VERSUS BAFFLE DEPTH FOR 
LIQUID AMPLITUDE g‘O.SOQ WITH THE BAFFLE 

WIDTH AS A PARAMETER (MILES FORMULA) 
REVISED 


* 



Figure 5- AO 
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DAMPING RSCTOR VERSUS BAFFLE DEPTH FOR WIDTH 
W» 0.025a WITH LIQUID AMPLITUDE AS A PARAMETER 
(MILES FORMULA) 


REVISED 






DAMPING FACTOR VERSUS BAFFLE DEPTH FOR WIDTH 
W >0.050 WITH LIQUID AMPLITUDE AS A PARAMETER 

(MILES FORMULA) 

REVISED 





Figure 5-42 
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DAMPING FACTOR VERSUS BAFFLE DEPTH FOR WIDTH 
WITM LIQUID AMPl/TUDE AS A PARAMETER 
(MILES FORMULA) 


REVISED 


v 

•$ 



Figure 5-43 
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DAMPING FACTOR VERSUS BAFFLE DEPTH FOR WGTH 
W-0.200 WITH LIQUID AMPLITUDE AS A PARAMETER 
(MILES FORMULA) 

REVISED 



Figure 5-44 
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DAMPING FACTOR VERSUS SURFACE AMPLITUDE FOR 

BAFFLE WIDTH W« 05a WITH BAFFLE LOCATION AS 
A PARAMETER (MILES FORMULA) 

REVISED 


rs 



Figure 5-45 
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DAMPING FACTOR VERSUS BAFFLE WIDTH FOR 

C w » IOa WITH LOCATION OF BAFFLE AS A PARAMETER 
(MILES FORMULA) 

REVISED 



Figure 5-47 
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DAMPING FACTOR VERSUS BAFFLE WIDTH FOR 


£ s ,30a WITH LOCATION OF BAFFLE ASA PARAMETER 
W (MILES FORMULA) 

REVISED 



Figure 5-48 
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DAMPING FACTOR VERSUS BAFFLE DEPTH TOR 
VARIOUS BAFFLE WIDTHS AND A LIQUID AMPLITUDE 



Figure 5-49 
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DAMPING FACTOR VERSUS BAFFLE DEPTH FOR 
VARIOUS BAFFLE WIDTHS AND A LIQUID AMPUTUDE 
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DAMPING FACTOR VERSUS BAFFLE DEPTH FOR VARIOUS 

BAFFLE WIDTHS AND A LIQUID AMPLITUDE OF £ »0.2a 

w 

AND A BAFFLE DISTANCE FROM EACH OTHER D-O.20 





Figure 5-52 


w 


.20 



.32 * * DAMD.MQ FACTOR VERSUS BAFFLE DEPTH FOR VARIOUS 

I BAFFLE WIDTHS AND A LIQUID AMPLITUDE OF C «0.50o 

j w 

I | AND A BAFFLE DISTANCE FROM EACH OTHER D»0.2o 



Figure 5-53 
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DAMPING FACTO** VERSUS BAFFLE DEPTH FOR VARIOUS 
LIQUID AMPLITUDE AND A BAFFLE WIDTH OF w«. 05 a 
AND A BAFFLE DISTANCE D«0.2o 



Figure 5-55 
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DAMPING FACTOR VERSUS BAFFLE DEPTH FOR VARIOUS 
LIQUID AMPLITUDE AND A BAFFLE WIDTH OF w.lOo 
AND A BAFFLE DISTANCE D-0.2o 




Figure 5-57 
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DAMPING FACTOR VERSUS SURFACE AMPLITUDE FOR BAFFLE WIDTH W=.05o 
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Figure 5-58 







Figure 5-59 









Figure 5-60 
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DAMPING FACTOR VERSUS BAFFLE WIDTH FOR C w "- 3 ° 
W’TM LOCATION OF BAFFLE AS A PARAMETER 



Figure 5-61 
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DECREASING LIQUID LEVEL 








BAFFLES ABOVE AMO 
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DECREASING LIQUID LEVEL 









where h = depth of liquid 
d = depth of baffle 

a = radius of right circular cylinder 
w — ring width 
c = radial clearance 

References : (35) 

Investigations by Silveira, Stephes and Leonard (35) show that 
for a given baffle depth: 

1. With the radial clearance constant: 

(a) The damping factor increases as the baffle width increases (Figure 
5-11) 

(b) The frequency increases as the baffle width decreases (Figure 5-16). 

2. With the baffle width constant: 

(a) The damping factor decreases as the radial clearance increases 
(Figure 5-11) 

(b) The frequency increases as the radial clearance decreases 
(Figure 5-16). 
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5. 1.3.4 Conic Section 



d = depth of conic section below undisturbed free- 
surface 


a = radius of right circular cylinder 
w = conic section width 

References ; (35, 42, 43, and 83) 

Investigations by Silveira, Stephens and Leonard (35) show that 

1. For a given baffle depth: 

(a) The damping factor of a perforated and nonperfor ated conic section 
increases as the baffle width increases (Figure 5-68). 

(b) The frequency for a perforated and nonperforated conic section 
increases as the baffle width decreases (Figure 5-64 and 5-67). 

2. Although the damping provided by the conic sections is slightly higher 
than that for the ring baffles, the surface area of a conic section 
having the same baffle width as that of the ring baffle is considerably 
higher. 


Investigations by Liu (83) give a conformal mapping solution 
for the pressure distribution on the conic section (article 5. 1.7. 3). 
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Variation of 


section 



.2 4 .6 


Baffle locafion, -^s 


r with baffle location for conic- 
baffle. (35) 


e 5-64 



Variation of frequency with baffle location for 
perforated-conic-section baffle. (35) 


Figure 5-65 
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Variation of damping factor with baffle location for conic- 

section baffle. R = 6 inches: — = 2. (35) 

R 

i 

Figure 5-66 
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Variation of damping factor with baffle location for 

perforated-conic-section baffle. R = 6 inches; — = 2. (35) 

R 


Figure 5-67 
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5 . 1 . 3 . 5 Inverted Conic Section 



a = radius of circular cylindrical tank 
d g = depth of baffle 
w = inverted conic section width 


References ; (35, 42, and 83) 

of Investigations by Silveira, Stephens and Leonard (35) show the 
effect /baffle width on both the frequency and the damping factor as a func- 
tion of baffle depth (Figures 5.68 and 5.69). 

Investigations by Liu (83) give a conformal mapping solution 
for the pressure distribution on the baffle. 
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-.2 -.1 0 .1 .2 .3 A .5 .6 .7 

Bcffle location, -j^-S 

Variation of frequency with baffle location for inverted- 
conic-section baffle. (35) 


Figure 5-68 
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Baffle location, 


Variation of damping factor with baffle location for 
inverted-conic-section baffle. R = 6 inches; = 2. (35) 


Figure 5-69 







5.1.4 Cruciform Baffle 




where h = depth of liquid 

a = radius of right circular cylinder 
w = cruciform width 


References ; (35) 

Cruciform baffles are most effective against rotary slosh which 
has large horizontal components of fluid velocity near the tank wall. Cruci- 
form baffles are also slightly effective against lateral slosh since large 
horizontal fluid velocities exist where the node line of the first slosh mode 
intersects the tank walls. 

Investigations by Silveira, Stephens and Leonard (35) indicate, 
as would be expected, that damping is independent of fluid depth except for 
shallow fluid depths. 
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5. 1.5.1 Mat 



References : (40, 42, 43) 


Investigations by Eulitz (40) have been made using commercial 
coco-fiber mats cut into a shape corresponding to the cross-sectional area 
of the tank. Several layers were used to attain the required thickness and 
hollow, airtight aluminum spheres were inserted between the layers to pro- 
vide the necessary buoyancy. The device was highly effective in damping the 
liquid amplitude; however, because of the following, it is impractical for 
missile application. The inner walls of a propellant tank are usually 
obstructed by stiffener rings, pipelines and other equipment. Therefore, 
if a floating device is to be used, it must be capable of adapting to the 
changing cross-sectional area of the container. 
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5. 1.5. 2 Can 


Liquid Surface 



•Cans (simplified) 


References : (40, 42, 43) 


Investigations by Eulitz (40) have been made using cans made 
of perforated sheet metal and enclosing a hollow air-tight metal ball of 
slightly smaller diameter, which provides the necessary buoyancy. The 
optimum dimensions of a can necessary for adequate buoyancy and thereby the 
total number of cans is determined in the following manner. The fraction 
of submerged volume of the enclosed ball A is given by 

X = < ^HE^XS)[3 + ( 3 + n ) P] 


where 6 
Y 


n 

D 

P 


the thickness of the can wall 

the specific gravity of the can material 

the number of cans along a tank diameter 

the tank diameter 

the percent perforation 


For a particular set of values 6, y and P, a graph of X vs. D can be plotted 
consisting of a family of curves corresponding to various values of n (for ! 
example, see Figure 5-73). The curve passing nearest to the intersection of 
the desired mean values of X = 0.5 and the tank diameter D under considera- 
tion will give the optimum value of n. The total number of cans is then 
determined by 

N = i (3n 2 + 1) 
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5. 1.5.2 Can (continued) 


and the can diameter is given by 



The cans are made of perforated aluminum sheet enclosing a 
hollow air-tight aluminum ball, of slightly smaller diameter, which provides 
the necessary buoyancy (Figures 5-70 and 5-71). The size and number of cans 
used in the tests are indicated by n = 7 or N = 37 (see Figure 5-72). 

The containers are rigid-wall circular cylindrical tanks of 
diameter 17.5 and 25 inches, respectively. 

The optimum dimensions of the can necessary for adequate 
buoyancy, and thereby the total number of cans, is determined as follows: 
the desired value of X the fraction of submerged volume of the aluminum 
ball, is 0.5. Figure 5-73 gives an equation for X as a function of the 
thickness of the can wall , the specific gravity of the can material Y, 
the percent perforation P, the tank diameter D and the number of cans along 
a tank diameter n. The curves shown in Figure 5-73 enable the designer to 
determine the values for n and N corresponding to ^ 2 0,5 as a function of 
D for a set of values of 6, Y and P. Note that the notation used in 
Figures 5-71, 5-72, and 5-73 is applicable only to these figures. 
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D • TANK DIAMETER 

d|« DIAMETER OF THE CYLINDRICAL BODY 


dt*DIAMETER OF THE CONE OPENING 
d a * DIAMETER OF THE ALUMINUM BALL* 
• * SLANT HEIGHT OF THE CONE 
J • WALL THICKNESS OF THE CAN 
% • SPECIFIC GRAVITY OF THE CAN 
MATERIAL 



0 * SURFACE OF THE CAN 
V * VOLUME QF THE CAN MATERIAL 
W* WEIGHT OF THE CAN MATERIAL *0*/ *y 
Ocome ■ ^d|*TT (eq.5) 

Ocyl « xdi 2 Trn (eq.4) 
n • number of bodies upon d 

0«ody*Ocyl + 20cone *4 TT(3+n)p 

(eq.6) 


P • PERFORATION IN % 

^SPHERE ■ -|-d| 2 TT (eq.7) 

^CAN * 0 BOOY + 0 SPHERE 

• 4*-7T [3 + (3+n)p] teq.8) 

W CAN * ^TT |5 i* (3+n)jj|*y (eq.9) 


Typical dimensions of cans. (40) 


Figure 5-71 



Tank Wall 




nd • 4d 


7d » D 

Total Number of cans 

Hexagonal Configuration of the cans | 
upon the surface. 


Liquid Surface' 



Figure 5-72 Distribution of cans within a tank. (40) 
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Figure 5-73 Optimum Number of Cans for Circular Cylindrical Tank. (40) 







5*1*6. 1 Cross Partitions 



References ; (4, 58) 

Investigations by Bauer (4) indicate that cross partitions have 
a greater effect on the eigen frequencies than do concentric partitions. For 
the quarter-sectored tank, the oscillating propellant takes on various modes 
of vibration. Also, in the case of a quarter-sectored tank, the vibrating 
liquid mass is reduced to more than one-half of that of a circular cylindri- 
cal tank. 


Investigations by Garza (58) indicate the effect of hole size 
on the frequency as a function of percent perforation (Figure 5-78 and 5-79) 
and the effect of both percent perforation and frequency on the damping ratio 
as a function of the excitation amplitude (Figure 5-80 and 5-81). 
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5. 1.6.2 Concentric Partitions 



References ; (4, 52) 



5.1.7 Pressure Distribution . The designer of a baffle system is 


interested in finding the pressure history of the fluid acting on the baffles 
because after the spacing and width of the baffles have been optimized, 
integration of the pressure will give the moments to be expected during 
flight. Thus, with the forces and moments acting on the baffle determined, 
the baffle material and its thickness can be chosen to withstand the load 
exerted on the baffles by the oscillating fluid, insuring, also, that the 
system pays a minimum weight penalty from the addition of the baffles. 


5. 1.7.1 Procedure I (28) . The force per unit area on the ring 


is given by 


p(0 ,t) = • y P w 2 


= C D Tj- 0 J 2 f 2 (-d) cos 2 0 sin 2 (o)t) 


(5-9) 


where w is the vertical component of the fluid velocity, f(z) = sinh k(z+h)y 
sinh(kh) and C D is the local drag coefficient. Its value may be obtained 
from the empirical relation 


C D " 15 XJ 


U T 
m 

D ’ 


U T 
m 


— D — 


< 20 


and 


C D = 2 ’ 


U T 

-g- > 100 


U T 

where is the "period-parameter"; U denotes the timewise maximum 

D m 

velocity, T the period, D the plate width, and 


U m T [° J S 1 f (- d ) cose ] (2ir/w) 2*f(-d) / 

— — = V — ) cos 0 


aa 


5-107 



2 

Substituting w = kg tanh (kh) into Eq. (5-9) and expressing the maximum 
pressure on the baffle as an equivalent head of the liquid, we have 


Z 

max 


P 

max 

Pg 


= ~ k a tanh (kh) f^(-d) (C ) 

2. u max 



( 5 - 10 ) 


5. 1.7.2 Procedure II (18) . Analytical research has 
shown that the motion of the free surface of tanked liquid propellants is 
normal to the undisturbed liquid free surface. Thus, the maximum pressure 
and moments exerted by the liquid on the baffles occur when the mean free 

surface is at baffle level. That is, when the baffle depth, with respect to 
the free surface, is zero (d = o). 


The dynamic pressure of the fluid acting at any point on a baffle 
in a cylindrical tank can be expressed as 


P 


2kgC ^ p cos 6 
144 (e^ - 1) 


V'l f> 


V'I> 



(5-11) 


where 


p = dynamic pressure, psi 

2 

kg = longitudinal acceleration, ft /sec 

3 

P = mass density of the liquid, #m/ft 
a = radius of the cylinder , ft 

= angle corresponding to sloshing amplitude 
^1 = a maximum amplitude of sloshing? ft 
0 = angular coordinate 
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5 • 1 • 7.2 Procedure II (18) (continued) . 


It can be seen from Equation (5-11) that the maximum positive pressure and 
maximum negative pressure will occur at 0 = o, r=a and 0 = 180°, r=a, 
respectively. The analytical expression for the moments acting on the baffle 
is obtained from the dynamic pressure. It has the following form: 


V s) = 


7 P (aV)S 2 X(S) - l 


sSxs y 

n=l V S z -k/ 


cosh (— )(h-d) 
a 

e^(e^-l )cosh(~^) 


n /L n n 


K(— ,e ) 
a n 


(5-12) 


where 


a n 


a 

i - (— ) 

a 


e a 

v-v> 


W , 


+ <— ) 2 <e ) 

a n 


. e a « 

W 


= constant 


S = Laplace transformation operation 

a = Inner radius of annular baffle, ft 
o 


For a circular cylindrical tank of arbitrary size it is convenient to express 
Equation (5-12) as a ratio of the baffle moment to the product of longitudinal 
acceleration kg and sloshing amplitude 5^. That is, 


I PL 


e i d 


pa k|c 1 


ge 


K(-£, e ) 

a l 1 


(5-13) 


where 


Im^ I = absolute value of the moment on the baffle corresponding to 
the first oscillation mode (n = 1) 
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5 . 1 . 7 . 2 Procedure II (18) (continued) 


Graphical representation of Equation (5-13) can be found in Figure 

(5- 74 ). It shows the baffle moment as a function of baffle width. The 

curves in this figure are valid only for fluid depth greater than one radius 

h h 

of the tank measured from the bottom. That is, — > 1* For — < 1, the 

a a 

following correction factor is recommended: 

V 

CF FJd e^d Ejd 

cosh + sinh tanh 

a a a 


Then, the corrected moment expression becomes 



(5-14) 
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Moments Produced on Ring Baffles (18) 


Figure 5-74 
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5. 1.7.3 Procedure III (83) . Conformal mapping is used 


to transform the geometrically complex representation of the baffle and the 
tank wall into a simple plane for which a solution to the Bernoulli equation 
for incompressible, nonviscous, unsteady flow is readily obtainable. Since 
this equation gives the velocity and pressure distribution on the baffle 
in the transformed plane it is then a simple matter to take the inverse trans- 
formation to find the solution in the original plane. 

The transformation is obtained by using the Christoffel- 

Schwarz expression defined as 

0 -1 o -l CL-1 

z ~ A j (w-u 2 ) ... (w-u k ) k dw + B (5-15) 

where 

z = F(w) = x + iy 
w = u + iv 

x,y = Cartesian coordinates in z-plane 

u,v = Cartesian coordinates in w-plane 

A = Constant which determines the scale factor and rotation of 
the polygon 

B = Constant which determines the origin of the w-plane 

Tra^ = Interior angles of the polygon formed by the tank wall and 
the baffles in the z-plane 

= Vertices of the polygon on the u-axis 

The pressure distribution applicable to all three cases 
herewith given is as follows: 
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5. 1.7.3 Procedure III (83) (continued) 


p ( u * v ) = P £ |u|ui [rpF(w) - Auj - |u| 2 (| - 2ARf + |^)} (5-16) 

where w is the circular frequency, P the density of the fluid, U the 
velocity of the fluid acting on the baffle and f, K, and R are parameters 
depending on the geometry of the system. 

Baffle Type I. The baffle-tank wall system in the z-plane is map- 
ped onto the w-plane by the transformation 

w 

2 = A f (w+c) (w-u 2 )" 1//n (w-u 3 )(w-u^) 1 ^ n (w-c)" 1 ^ 2 dw + B (5-17) 

' o 

The result is shown in the figure for Baffle Type I below. 

i 

i 


z-plane w-plane 

y 



Baffle Type I 
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5. 1*7.3 Procedure III (83) (continued) 


The parameters and constants for the pressure distribution 
given in Equation (5-16) are as follows: 


K(u,o ' 


j/lk 


2(c 2 -u 2 ) 


u-u 


u-u 


4 

2 


1/n 


R(u,o) 


2A(u-u 3 ) 2 


u-u 


u-u 


2 

4 


1/n 


f(u,o) = 8 sin(6 3 /n) 


g(ujo) — -8 cos ( © 3 /n ) 

A = b/I 2 B = + ib^ 


0 3 (u,o) 


o for -c < u < u 2 » < u < c 

TT for U 2 £ u < u^ 


8(u,o) 


r. 


2(e S 77 


<u-u 3 ) 


I = 
o 



o 


|F' (w) |dw 


3 

I 2 = f (w+c) 1 ^ 2 (w-u 2 > (u 3 -w)(u. -w) 1 ^ 11 (c-w)" 1 ^ 2 dw 

U 2 

b^ = -b(l+l o /l 2 > sin a ; b 2 = b(l+I Q /I 2 ) cos a+ a 
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The constants u^> can be determined from the 

integral expression for I above and from the following definite integrals 

I = f (w+c)' 1/2 (u 2 -w)" 1/n (u 3 -w)(u 4 -w) 1/n (c-w) 1/2 dw 
* -c 

1 3 = f (w+c)" 1 ^ 2 (v-u 2 ) 1 ^ n (w-u 3 )(u 4 -w) ^ (c-w) ^ dw 

u 3 (5-18) 

1 4 = f 5 (w+c)‘ 1/2 (w-u 2 )" 1/n (w-u 3 )(w-u 4 ) 1/n (c-w )" 1/2 dw 

U 4 


where 1^, I^* I 3 , and I 4 are definite integrals which are proportional to 
the length of the sides of the polygon ABCDE corresponding to polygon 
12345 in the transformation for Baffle Type I. Thus, 



The unknown constants u 3 , u 4 can be determined by finite difference 
techniques by using the simultaneous equations 


I k (u 2 , u 3 , u 4 ) = X i c I i^ u 2 ,u 3 ,u 4^ » k = 2, 3, 4 


(5-19) 


For example: let u^ 9 u^ and be the solution of the above equation and 

take u^°\ u^°^ and uf°^ as the initial values differing little from the 


exact solution u^, and u^. We expand equation (5-19) in a Taylor f s 

(1) * (o) 

~ u k " u k anC * terml - nate at t ^ ie 


series in powers of the differences 
first order terms, thus obtaining the first approximation 
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5. 1.7.3 Procedure III (83) (continued) 



+ I 

j — 2 , 3 ,4 





+ l 

j=2,3,4 


a j (1)l l 


(o) 1 

1>U j] 


k— 2*3,4 (5_20 ) 


where 



The partial derivatives which have limits not depending on u 2> u 3 and 
may be obtained by direct differentiation, i.e., performing the partial 
differentiation, then integrating. Those that have u 2> u 3 and u 4 as 
limits are computed numerically by using 


T (m) 

k -“j 


( 


V u 


(m) 

j 


+ Au) 



(5-21) 


where Au is a small increment. Solving for simultaneously from 

equation (5-20) and using the second approximation, 


u 


( 1 ) 

j 



+ a 


( 1 ) 

j ’ 


j=2, 3, 4 


in equation (5-20), we obtain 



+ l 

j=2,3,4 






+ l 

j=2,3,4 




k=2 ,3,4. 


The above process is carried out until 



s converge. 
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5. 1.7.3 Procedure III (83) (continued) 


Special Case: Baffle Type IA . This is a special case 

of Baffle Type I for which a = 0, u 2 = -c, and u^ = c. The transformation 
becomes 

w _ ( n+ 2) _ (n-2) 

z = A f (c+w) (w-u_ ) (w-c ) dw + B (5-22) 

' o 

The results are shown in the figure for Special Case la below. 



For this case, the parameters and constants needed to obtain 

the pressure distribution (Equation 5-16) are the following: 

_ (n+1 ) _ ( n -l ) 

K(u,o) = A (c+u) n (c-u)’ n 


R(u,o) = 


/2 

2A 


c+u 

C-u 


1 1/n 


(u-u 3 )' 


f(u,o) - - / 2(c^-u^) (u-u 3 )sin o 


(5-23) 


g(u,o) - - /^7 ^u-u 3 )cos o 
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5. 1.7. 3 Procedure III (83) (continued) 


A = b/I 2 
B = + ib 2 

ReF(u,0) = Asina [ (c + u) ^ n+2 ^ 2n ( u 


V 


(c - u) 


>(n-2)/2n 


du + b, 


ImF(v,0) = - Acos a 


f (c + u )-^ n+2 >/ 2 " 


(u - u^) (c - u) 2 ^ 2n <j u 4 b 2 


The pressure distribution for Case IA is shown in figure (5-75). 

Special Case; Baffle Type IB . This is a special case of Baffle 
Type I for which ot = o, u^ = -c, u^ = o, u^ = c, and ^ = o. The transformation 
becomes 


z = A 



wdw 



= A 



(5-24) 


The results are shown in the figure for special Case IB below. 


y 

3 


tank wall "77777/7 7TTT7777 




5 


c u 

T~ 


Special Case lb 

The pressure distribution is well known. It is parabolic in shape v ..ch the 
maximum disturbance pressure occurring at the base of the baffle and the 
minimum (z ero) at the tip of the baffle. 
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5. 1.7. 3 Procedure III (83) (continued) 


Baffle Type II . The baffle-tank system (n^ = = n, 

u 9 = o) in the z-plane is mapped onto the w-plane by the transformation 


z = A 


r 


(c+w) 


<n r 2) 

2 n, 


1 - — 
(w-u 2 > ni 


(n 2 - 2 ) 


l 2 (w-c) 


2 n 


2 dw + B (5-25) 


The results of this transformation is shown in the figure below. 

w-plane 



rank wall 


3 x 

/ / / / rr 


AU 

« 

-c 


u a o 


1 


c 


Baffle Type II 

The parameters and constants given in Equation (5-16) can be obtained as 
follows : 


P- = AK- (sin a. + icos a, ) = «£ + ig.) 
dw ii 

^ = —77 (sin a. - icos a, ) = R(f - ig) (5-26) 

dz AK 1 1 i. 

- c < u < u 2 
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5. 1.7. 3 Procedure III (83) (continued) 


— = AK' (sin a 2 - icos a 2 > = K(f + ig) 


^ (sin o 2 + icos a 2 > = R(f - ig) 

” u 2 1 u i c 


(n 2 -2) 

K' = (c+u) 2n 2 


u-u. 



(n 2 -2) 
(c-u) 2n 2 


A = 


b 

I^cosa^ 



K'du 


b^ = b(tan a 2 


sin o„ I 

2 _£ ) 

cos 1^ 


I 

c 


i 


K f du 



cos a 2 
cos ct^ 



The constant can be found from 


cos a. 

I, * I - I 0 , 

cos a 2 1 c 2 




K'du 


(5-27) 


(5-28) 


5-120 


5.1 # 7. 3 Procedure III (83) (continued) 


Special Case: Baffle Type IIA . Same as Baffle Type II 

but ot = ot^ = ^2* T ^ ie trans ^ ormat ^ on becomes 

w _ (n-2 ) (n-2) 

z = A | [(w^-c^) (w) n ] dw + ib (5-29) 

' o 


The results are shown in the figure below 


2 



Baffle Type Ila 


Baffle Type III . The transformation from the z-plane 
to the w-plane is given by 


z 


= A 


i 


(w 2 -^) dw 
} (w 2 -c 2 )(w 2 -u 2 )' 


(5-30) 
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5. 1.7. 3 Procedure III (83)(contlnued) 


and the results shown in the figure below. 
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BAFFLE TYPE III 

The parameters and constants given in Equation (5-16) are as follows: 


— = a 

dw 


, 2 2 v 

(u -u 2 ; 


dw __ _ 
dz A 


[71 2771 2 ? 

y(c -u )(u -u^) 

L 2 2 W 2 2,\ 

y Ce -u Xu -u 3 ; 


(o - i) = K(f + ig) 


2 2 (o + i) = R(f - ig) 

(u -u 2 ) 


(5-31) 


for u^ < | u | < c 


dz = A 

dw 


( 2 2 \ 
(u -u 2 ) 


dw _ _ 
dz A 


J, 2 2 W 2 2 ' 

l(c -U Xu^-u ) 

[ 71 - 1777-27 

J ^(c -U Xu^-U ) 


2 2 
u -u 2 


(1 + io) = K(f + ig) 


(1 - io) = R(f - ig) 


(5-32) 


for | u | < u„ 
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5. 1.7. 3 Procedure III (83) (continued) 


The constants u^, u^ » 


from the equations 


Kl /K 2 = d/2b , K, = K 3> 
r U 3 (u 2 -u 2 ) du 

K l = J 


«2 = [ 


[71 2w 2 2 ? 

'j (c -u Ku^-u ) 

u 2 (u 2 -u 2 )du 

u„ J. 2 2w 2 2? 

3 ](c -u )(u.,-u / 


K 3 = 


-j: 


(u 2 -u 2 )du 


(c 2 -u 2 )(u 2 -u 2 ) 


\ 


and A are determined simultaneously 


A = d/2K 1 , B = o 


(5-33) 


Special Case: Baffle Type IIIA . The results are ident- 

ical with Special Case Baffle Type IB but with u^ = u 2 = 0. 
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INCHES 


U 



PRESSURE ON A SUNT BAFFLE (83). 


Figure 5-75 
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5-1.8 Conclusions and Recommendations 


In choosing a damping system, the designer should start with 
solid ring baffles and use Miles ! expression, Equation (5-1), for his 
first approximation of the damping factor. The reason is that Miles 1 
equation is simple in format and in nature and a rather rapid "ball-park 
figure” for the spacing and width of the baffles can be obtained. Bauer's 
extension (article 5. 1.3. 2. 3) would then give more accurate results and 
optimum conditions. The pressure acting on the baffle (Equations 5-10, 5-11, 
and 5-16) will permit the optimization of the thickness and material of each 
baffle and thus the calculation of the total weight added by the damping 
system. 

Results 

1. The maximum frequency occurs when the baffle is located at the 

d s 

undisturbed free liquid surface, i.e., — = 0 (Figures 5-2 and 5-3). 

Note: The curves in Figures 5-2 and 5-3, for high values of 

d s 

— , approach asymptotically, the first resonant frequency for a 

circular cylindrical tank containing no baffles. 

2. Agreement between experimental values of the damping factor and 

those obtained from Miles' equation is good except in the range of 
d g 

0 < ^ < 0.125 (Figure 5-5). 

3. With the hole size constant and for a given baffle depth: 

(a) The frequency increases for increasing percent perforation 
(Figure 5-2). 
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Results (continued) 


(b) The damping factor increases for decreasing hole size 
(Figure 5-8). 

4. With the percent perforation constant and for a given baffle depth 

(a) The frequency increases for increasing hole size 
(Figure 5-3). 

(b) The damping factor increases for decreasing hole size 
(Figure 5-7). 

5. The frequency of a solid fixed ring baffle is higher than that of 

d s d s 

a perforated baffle for — $ 0.05 and lower for * 0.05 

(Figures 5-2 and 5-3). 

6. Damping produced by perforating baffles is consistently lower 

than that produced by a solid baffle (Figure 5-7 and 5-8). 

7. The damping factor increases as the excitation amplitude increases 

(Figure 5-6). 

8. For a given baffle depth: 

(a) The damping factor of a fixed ring, conic section and per- 
forated conic section Increases as the baffle width increases 
(Figures 5-10, 5-67, and 5-68). 

(b) With the radial clearance constant, the damping factor of a 
ring-with-radial-clearance increases as its baffle width 
increases (Figure 5-11). 

(c) With the baffle width constant, the damping factor of a ring- 

* 

with-radial-clearance decreases as the radial clearance 
increases (Figure 5-11). 
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9. The damping factor is independent of kinematic viscosity and 

increases as the excitation amplitude increases (Figures 5-13 and 
5-14). 

10. Agreement between experimental values of the damping factor and 

those obtained in Miles f equation is good except for 
de 

0 < ^ < (0.8 - 1.2), (Figures 5-12 and 5-14). 

11. For a given baffle depth: 

(a) The frequency for a fixed ring, conic section, and perforated 
conic section increases as the baffle width decreases 
(Figures 5-15, 5-64, and 5-66). 

(b) With the radial clearance constant, the frequency increases 
as the baffle width decreases (Figure 5-16). 

(c) With the baffle width constant, the frequency increases as 
the radial clearance decreases (Figure 5-16). 

12. During propellant loading, the cans erected smoothly and arranged 
themselves in the required configuration even though they were 
randomly placed in the tank bottom before starting the filling 
operation. For this reason, they operated as good anti-vortex 
devices . 

13. The devices floated perfectly and yielded high damping effects. 

14. Draining the liquid while the tank was oscillating did not change 
the damping effect. 

15. An inherent disadvantage of cans is that they are free to move when 
the container is empty, i.e., they are loose bodies during the 
transportation of the missile and can thereby cause damage to the 

laasr tank 
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16 . 


The dashed parts of the curves in Figures 5-9, 5-10, and 5-67 
indicate that, because of the turbulent nature of the liquid in the 
regions of higher damping, an uncertainty existed in defining the 
exact location of maximum damping. 

17. A divided baffle which extended down to the ring baffle was used 
for most of the tests in Reference (35). This tended to increase 
the value of the damping as seen in Figure 5-14. 

Conclusions 

1. Fixed ring baffles can be perforated with holes of relatively small 
diameter, thereby reducing the baffle area by as much as 23% with 
no appreciable loss in damping effectiveness. Such reduction in 
baffle area may be effectively used for additional baffles thus 
increasing the minimum damping ratio value with no appreciable 
increase in baffle weight. 

x o 

2. For an excitation amplitude — = 0.00417, the baffle above the 

d s 

liquid surface adds to damping only when the distance — < 0.125 

and the baffle below the surface contributes to damping only when 
d s 

0 < — < 0.375. The baffle spacing for some minimum Y can be 

— K ~ s 

determined by adding the depth of the submerged baffle to the depth 
d s 

of — = 0.125 for which the upper baffle is still effective. 

3. Fixed ring baffles can efficiently introduce liquid damping into 
propellant tanks. 

4. Fixed ring baffles are not very efficient for shifting the liquid 
resonant frequency. 
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Conclusions (continued) 


5* Based on the total surface area of a baffle, the highest mean 

damping factor of the baffles tested appears to be afforded by the 

fixed ring baffle (Figure 5-78). 

Summarizing, the following remarks are pertinent. 

1. Bauer f s extension to Miles f equation can be used without 
further experimental verification except for a check of the final product. 

2. Ring baffles provide little frequency effect. 

3. Solid ring baffles with small diameter perforations (about 20% 
open area) can effectively be used to reduce baffle weight addition without 
an appreciable loss in damping effectiveness. 

4. Conical rings offer greater degree of damping than floating 
cans, except where the cans carefully cover the liquid surface. 

5. The amount of damping produced by both floating cans and conical 
rings is strongly dependent upon excitation amplitude. 

6. Cans are the only suitable type of float. However, ring baffles 
are more reliable and appear to give greater damping. 

7. Floats, like baffles, are used mostly to damp liquid motion 
however, floats are like partitions, in that their effect on frequency and 
damping is constant for all shallow depths. 

8. Subdivision by radial walls is much more effective in reducing 
the amount of sloshing than subdivision by concentric walls. In the case 

of a quarter-sectored tank the vibrating liquid mass is reduced to more than 
1/2 that of a circular cylinder. 

9. Ring baffles are far better as damping devices than conic, 
inverted conic or cruciform baffles. 
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10. Cross partitions have a greater effect on the eigen-frequencies 
than do concentric partitions. 

11. Floats should cover the entire cross-sectional area and fill 
the zone of oscillation, i.e., to a depth of one-fourth the tank diameter. 

12. Floats should be capable of following the changing liquid level. 

13. Floats should be capable of adapting to the changing cross- 

* 

sectional area of the container, i.e, the float should not catch-on or cling 
to stiffener rings, pipelines or other equipment obstructing the inner walls 
of the tank. 

14. Floats should not change the moment of inertia of the liquid 
when subjected to roll excitations about the longitudinal axis of the con- 
tainer. 

15. Partitions should divide the total liquid mass into approximately 
equal partial masses thereby effectively reducing the tank diameter. 
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(d) Inverted conic 
section. 


(e) Perforated 
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Baffle configurations. (35) 


Figure 5-76 
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Effect of hole size and percent perforation 
on resonant frequency for a quarter sector tank. (58) 


Figure 5-78 
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Effect of hole size and percent perforation 
on resonant frequency for an eighth sector tank. (58) 


Figure 5-79 
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Effect of percent perforation and resonant 
frequency on the damping ratio as a function of 
excitation amplitude for a quarter sector tank. (58) 


Figure 5-80 
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TABLE A- 3 


ZEROS r (m=l) OF THE FIRST DERIVATIVE OF THE BESSEL 
ran 

FUNCTION OF THE FIRST ORDER AND FIRST KIND = OJ. 


n 

«n 

n 

«n 

n 


0 

1.84119 

a 

36.89000 

22 

71.45677 

1 

5.33144 

12 

40.03341 

23 

74.60000 

2 

8.53631 

13 

43.17665 

24 

77.74383 

3 

11.70600 

14 

46.31957 

25 

80.88514 

4 

14.86359 

15 

49.46238 

26 

84.02718 

5 

18.01553 

16 

52.60507 

27 

87.16914 

6 

21.16440 

17 

55.74758 

28 

90.31228 

7 

24.31145 

18 

58.89000 

29 

93.45301 

8 

27.45751 

19 

62.03234 

30 

96.59499 

9 

30.60194 

20 

65.17465 

31 

99.73675 

10 

33.74397 

21 

68.31682 
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TABLE A-5 

ROOTS c OF THB FIRST DERIVATIVE OF THE BESSEL FUNCTION 
urn 

OF ORDER 4m AND OF THE FIRST KIND [ J» (e ) = 0 3. 

Hin inn 



1 

3 

5 

0 

5.31756 

13.82109 

22.26759 

1 

9.28240 

18.74485 

27.71172 

2 

12.68191 

22.62927 

31.97366 

3 

■ 

15.96411 

26.24604 

35.87393 

4 

19.19603 

29.72898 

39.58453 

5 

22.40103 

33.13145 

43.17654 


A-5 




ROOTS c (m=l) OF THE DETERMINANT 



A -6 





















Roll about the z-axis 



A -9 




REFERENCES 


1. Chu, W-H., "Free Surface Condition for Sloshing Resulting from 

Pitching and Some Corrections", ARS Journal , (Nov. 1960), 
1093-1094. 

2. Schmitt, A. F. , "Forced Oscillations of a Fluid in a Cylindrical 

Tank Oscillating in a Carried Acceleration Field — A Correction", 
Convair Rep. ZU-7-074 , (Feb. 1957). 

3. Bauer, H. F. , "Theory of Fluid Oscillations in Partially Filled 

Cylindrical Containers", MSFC, Rept. No. MTP-AER0-62-1 , 

(Jan. 1962). 

4. "Fluid Oscillations in the Containers of a Space Vehicle 

and their Influence Upon Stability", NASA-TR R-187 . 

5. Chu, W. , "Sloshing of Liquids in Cylindrical Tanks of Elliptic Cross- 

Section", Southwest Research Institute , Technical Report No. 2 , 
Contract DA-23-Q72-ORD-1251 , (Sept. 1958). 

6. Okhotsimskii, D. E., "Theory of the Motion of a Body with Cavities 

Partly Filled with a Liquid", NASA-TT , No. F-33 . (May 1960)-- 
Translated from Prikladnaia Matematika : Makhanika , Vol. 20, No. 1 
(Jan. -Feb. 1956). 

7. Bauer, H. F. , "The Moment of Inertia of a Liquid in a Circular 

Cylindrical Tank", ABMA Report No, DA-TR-5-58 , (April 23, 1958). 

8. Zhukovskii, N. E., The Motion of a Rigid Body with Cavities Filled 

with a Homogeneous Flowing Liquid . Collected Works (in Russian) 
(Gostekhizdat, M.-L., 1949), Vol. II. Translated by Consultants 
Custom Translations, Inc., 227 West 17th Street, New York 11, N.Y 

9. Lamb. H. , Hydro dynamic s . Sixth ed. , Dover Publications, 1945. 

10. Rayleigh, Lord, "On Waves", Phil. Mag. , Series 5 , Vol. 1 , (1876) 

Seite 257-279. 

11. Bauer, H. F. , "Fluid Oscillations in a Circular Cylindrical Tank", 

ABMA Report No. DA-TR-1-58 , (April 18, 1958). 

12. Eulitz, W. , "The Sloshing Phenomenon and the Mechanism of a Liquid 

in Motion in an Oscillating Missile Container", Rep. No. DS-R-31 , 
POD , ABMA , (Oct. 1957). 

13. Abramson, H. N. , and Ransleben, G. E., Jr., "Simulation of Fuel 

Sloshing Characteristics in Missile Tanks by Use of Small Models" 
Tech. Rep. No. 3, Contract DA-23-072-ORD-1251 , SWRI , (March 1959) 

14. ., "Some Comparisons of Sloshing Behavior in Cylindrical Tanks 

with Flat and Conical Bottoms", Tech. Rep. No. 4 , Contract 
No. DA-23-072-ORD-1251 , SWRI , (May 1959). 


B-l 



15. "A Note on Wall Pressure Distributions During Sloshing in 

Rigid Tanks", Tech. Rep. No. 5, Contract No. DA-23-072-ORD-1251 , 
SWRI , (June 1959). 

16. McCarty, J. L. and Stephens, D. G. , "Investigation of the Natural 

Frequencies of Fluids in Spherical and Cylindrical Tanks", 

Langley Research Center , NASA-TN , D-252 , (May 1960). 

17. Abramson, H. N. , Chu, W. H., Ransleben, G. E., Jr., and Squire, W. , 

"A Theoretical and Experimental Study of Fuel Sloshing", Quarterly 
Progress Report No. 8 , Contract No. DA-23-Q72-ORD-1251 , SWRI , 

(May 1960). 

18. Koelle, H. H., Editor-In-Chief, Handbook of Astronautical Engineering . 

McGraw-Hill Book Co., Inc., New York, N.Y. , 1961. 

19. Bauer, H. F. , "Theoretical Investigation of Propellant Sloshing in a 

Circular Cylindrical Ring Tank", ABMA , Quarterly Report No. 1 , 
Redstone Arsenal, Alabama , (Dec. 20, 1959). 

20. Lorell, J. , "Forces Produced by Fuel Oscillations", Progress Report 

No. 20-149 , Jet Propulsion Laboratory , (October 16, 1951). 

21. Abramson, N. H. , and Ransleben, G. E., Jr., "Liquid Sloshing in Rigid 

Cylindrical Tanks Undergoing Pitching Motions", SWRI , Technical 
Report No. 11 , Contract No. DA-23-072-0RD-1251 , (May 1, 1961). 

22. Hutton, R. E., "Fluid Particle Motion During Rotary Sloshing", 

Space Technology Laboratories, Inc. , (March 1963). 

23. "An Investigation of Resonant, Nonlinear Nonplanar Free-Surface 
Oscillation of a Fluid", Ph.D. Dissertation, University of 
California , Los Angeles, California , (October 1962). 

24. Miles, J. W. , "Free Surface Oscillations in a Rotating Liquid", The 

Physics of Fluids , Vol. 2 (1959), 297-305. 

25. Miles, J. W. and Troesch, B. A., "Surface Oscillations of a Rotating 

Liquid", Journal of Applied Mechanics , (December 1961). 

26. Bauer, H. F. , "Fluid Oscillations of a Circular Cylindrical Tank 

Performing Lissajons-Oscillations", ABMA-TR , Huntsville, Alabama , 
(April 18, 1958). 

27. ., "Fluid Oscillations in a Cylindrical Tank with Damping", 

ABMA-TR , Huntsville, Alabama , (April 23, 1958). 

28. Miles, J. W. , "Ring Damping of Free Surface Oscillations in a Circular 

Tank", Journal of Applied Mechanics , (June 1958). 

29. Bauer, H. F. , "Tables and Graphs of Zeros of Cross Product Bessel 

Functions", MSFC Report No. MTP-AERQ-63-50 , (June 24, 1963). 


B-2 



30. Langner, C. G. , "A Preliminary Analysis for Optimum Design of Ring 

and Partition Antislosh Baffles", SWRI , Tech. Rept. No. 7 , 

Contract No. NAS8-1555 , Project No. 6-1072-2 , (April 30, 1963). 

31. Bauer, H. F. , "The Damping Factor Provided by Flat Annular Ring 

Baffles for Free Fluid Surface Oscillations", MSFC , Rept. No . 

MTP- AERO- 6 2-81 , (Nov. 13, 1962). 

32. Garza, L. R. and Abramson, N. H., "Measurements of Liquid Damping 

Provided by Ring Baffles in Cylindrical Tanks", SWRI , Tech. 

Rept. No. 5, Contract No. NAS8-1555, Project No. 6-1072-2 , 

(April 1, 1963). 

33. O'Neill, J. P. , "Semiannual Report on Experimental Investigation of 

Sloshing" , TR-59-0000-00713 , Space Tech. Labs., Inc. .(Jan. 1- 
June 30, 1959). 

34. Abramson, H. N. , Garza, L. R. , and Kana, D. D. , "Fuel Sloshing Studies", 

SWRI, Quarterly Progress Rept. No. 3 , Contract No. NAS8-1555 , 

(Jan. 1, 1962). 

35. Silveira, M. A., Stephens, D. G., and Leonard, H. W. , "An Experimental 

Investigation of the Damping of Liquid Oscillations in Cylindrical 
Tanks with Various Baffles", NASA TND-715 , Langley Research Center , 
(May 1961). 

36. Abramson, H. N. and Ransleben, G. E., Jr., "Some Studies of a Floating 

Lid Type Device for Suppression of Liquid Sloshing in Rigid 
Cylindrical Tanks", SWRI , Tech. Rept. No. 10 , Contract No. 
DA-23-Q72-ORD-1251 . 

37. "A Note on the Effectiveness of Two Types of Slosh Suppression 
Devices", SWRI, Tech. Rept. No. 6 , Contract No. DA-23-Q72-ORD-1251 , 
(June 15, 1959). 

38. Garza, L. R. , "A Brief Comparison of Ring and Asymmetrical Baffle 

Characteristics", SWRI, Tech. Rept. No. 6, Contract No. NAS8-1555 , 
(April 15, 1963). 

39. Cole, H. A., Jr., and Gambucci, B. J. , "Tests of an Asymmetrical Baffle 

for Fuel-Sloshing Suppression", NASA TN D-1036 , (July 1961). 

40. Eulitz, Dr. W. , "A Can-Type Anti-Slosh Device Derived from Basic 

Slosh Studies", ABM A , D.O.D. Rept. No. D SD-TR-4-58 , Redstone 
Arsenal, Alabama , (August 12, 1958). 

41. Cole, H. A., Jr., and Gambucci, B. J., "Measured Two-Dimensional 

Damping Effectiveness of Fuel-Sloshing Baffles Applied to Ring 
Baffles in Cylindrical Tanks", NASA-TN-D-694 , Ames Research 
Center, Moffett Field, Calif., (Feb. 1961). 


B-3 


42. Eulitz, W. R. , and Beduerftig, H. , "Experimental Investigation of 

the Responses of a Liquid in an Oscillating Container", Office 
of Ordnance Research, Proceedings of the Third Conference on 
the Design of Experiments in Army Research Development and 
Testing, Duke Station, Durham, North Carolina , (1958). 

43. Kachigan, K. , "Forced Oscillations of a Fluid in a Cylindrical 

Tank", Convair Rept. ZU-7-046 , (Oct. 1955). 

44. Schmitt, A. F., "Forced Oscillations of a Fluid in a Cylindrical 

Tank Undergoing Translation and Rotation", Convair Rept . 

ZU-7-069 , (Oct. 1956). 

45. Miles, J. W. , "On the Sloshing of Liquid in a Cylindrical Tank", 

Ramo-Wooldridge Rept. AM6-5 , (April 1956). 

46. Abramson, N. H., Chu, W-H., and Ransleben, G. E., Jr., "Representation 

of Fuel Sloshing in Cylindrical Tanks by an Equivalent Mechanical 
Model", SWRI , Tech. Rept. No. 8 , Contract No. DA-23-Q72-ORD-1251 , 
(June 30, 1960). 

47. Abramson, H. N. , and Ransleben, G. E., Jr., "Simulation of Fuel Slosh- 

ing Characteristics in Missile Tanks by Use of Models", SWRI, 

Tech. Rept. No. 7 , Contract No, DA-23-Q72-ORD-1251 , (April 25, 1960). 

48. Buckingham, E., "Model Experiments and the Form of Empirical Equations", 

Trans. ASME, Vol. 37 , (1915), 263-296. 

49. Bauer, H. E., "Mechanical Model of Fluid Oscillations in Cylindrical 

Containers and Introduction of Damping", MSFC, MTP-AERO-62-16 , 
(February 26, 1962). 

50. Abramson, H. N. , Martin, R. J., and Ransleben, G. E., Jr., "Application 

of Similitude Theory to the Problem of Fuel Sloshing in Rigid 
Tanks", SWRI , Tech. Rept. No. 1 , Contract No. DA-23-072-QRD-1251 , 

(May 23, 1958). 

51. Scanlan, R. H. , and Mendelson, A., "Structural Damping", AIAA Journal , 

Vol. 1, No. 4 , (April 1963). 

52. Bauer, H. E., "Mechanical Analogy of Fluid Oscillations in Cylindrical 

Tanks with Circular and Annular Cross Sections", MSFC, MTP AER0-61-4 , 
Huntsville, Alabama , (January 12, 1961). 

53. Miles, J. W. , "Stability of Forced Oscillations of a Spherical Pendulum", 

Quarterly of Applied Mathematics , (1962); also. Space Tech. Labs. , 

No. EM 11-2 , (March 1960). 

54. Stephens, D. G., and Leonard, H. W., "The Coupled Dynamic Response of 

a Tank Partially Filled with a Liquid and Undergoing Free and 
Forced Planar Oscillations", NASA TN, D-1945, Langley Research 
Center , (August 1963). 


B-4 



55. Abramson, H. N. , Garza, L. R. , and Kana, D. D. , "Fuel Sloshing 

Studies", SWRI , Quarterly Progress Report No. 2 , Contract 
No. NAS8-1555 , (Oct. 1, 1961). 

56. — - — ., "Fuel Sloshing Studies", SWRI , Quarterly Progress Rept. 

No. 1, Contract No. NAS8-1555 , (July 1, 1961). 

57. Abramson, H. N., and Ranselben, G. E. , Jr., "Total Force Response 

Resulting from Liquid Sloshing in a Rigid Cylindrical Tank 
with a Vertical Center Wall Baffle", SWRI , Tech. Rept. No. 9 , 
Contract No. DA-23-Q72-ORD-1251 , (May 1, 1961). 

58. Garza, L. R. , "Measurements of Liquid Natural Frequencies and Damping 

in Compartmented Cylindrical Tanks" , SWRI , Tech. Rept. No. 8 , 
Contract No. NAS8-1555 , (Feb. 5, 1964). 

59. Bowman, F. , Introduction to Bessel Functions . New York: Dover 

Publications. 

60. Walls, J. , Ryan, J. , and Mooney, J., "Torsional Pendulum Analogy 

for Quarter and Eighth Sectored Cylindrical Tanks Subjected 
to Roll Oscillation", Northrop Space Labs. , Tech. Memo. No. 9 , 

(Feb. 21, 1964). 

61. Abramson, H. N. , Garza, L. R. , and Kana, D. D., "Some Notes on 

Liquid Sloshing in Compartmented Cylindrical Tanks", SWRI , 

Tech. Rept. No. 1 , Contract No. NAS8-1555 , (Feb. 15, 1962); 
also, ARS Journal , (June 1962). 

62. Abramson, H. N. , Chu, W. , and Garza, L. R. , "Liquid Sloshing in 45° 

Sector Compartmented Cylindrical Tanks", SWRI , Tech. Rept. No. 3 , 
(Nov. 1, 1962). 

63. Bauer, H. F. , "Theory of the Fluid Oscillations in a Circular Cylindrical 

Ring Tank Partially Filled with Liquid", NASA-TN-D557 , (Dec. 1960). 

64. Mooney, J., Ryan, J., and Walls, J. , "Slosh Handbook", NSL Tech. 

Memo. No. 10 , (Feb. 21, 1964). 

65. Mooney, J., Walls, J., and Ryan, J. , "Liquid Natural Frequencies of 

a Tank with Partially Submerged Sector Walls", NSL Tech. Memo . 

No. 14 , (April 7, 1964). 

66. Bauer, H. F., "The Effective Moment of Inertia in Roll of Propellant 

and Roll Damping", ABMA Report DA-TM-67-59 . 

67. ., "Propellant Sloshing", Rept. No. DA-TR-18-58, C-AB-71091 , 

Dev. Operations Div. , Army Ballistic Missile Agency , Redstone 
Arsenal, Alabama , (Nov. 5, 1958). 

68. Ransleben, G. E., Jr., "A Theoretical and Experimental Study of Fuel 

Sloshing", Progress Rept. No. 1 , Contract DA-23-Q72-ORD-1251 , 
Southwest Research Institute, (May 1, 1958). 


B-5 


69. Keulegan, G. H. , and Carpenter, L. H., "Forces on Cylinders and 

Plates in an Oscillating Fluid", National Bureau of Standards, 
Report No. 4821 , (Sept. 5, 1956). 

70. Brown, J., "Laboratory Test of Fuel Sloshing", Rept. No. SM-Dev. 783 , 

Douglas Aircraft Co., Inc . 

71. Case, K. and Parkinson, W. C., "The Damping of a Liquid in a Right 

Circular Cylindrical Tank", Ramo-Wooldridge Rept. GM 45-75 , 

(Sept. 1956). 

72. Howell, E. , and Ehler, F. G., "Experimental Investigation of the 

Influence of Mechanical Baffles on the Fundamental Sloshing 
Mode of Water in a Cylindrical Tank", Ramo-Wooldridge Rept. 

GM 45, 3-87, GM-TR-69 , (July 1956). 

73. Langhaar, H. S., Dimensional Analysis and Theory of Models . John 

Wiley, 1953. 

74. Epperson, T. B. , Pengelley, C. D., Ransleben, G. E., Jr., Wilson, L. E., 

and Younger, D. G. , "Nonstationary Airload Distributions on a v 
Straight, Flexible Wing Oscillating in a Subsonic Wing Stream", 
(Title Unclassified), WADC TR-55-323 . (January 1956). 

75. Mooney, J., Basurto, E., Kincaid, J., Ryan, J., and Walls, J. , 

"Slosh Handbook", Progress Report No. 2, Northrop Space Labs. , 

Tech. Memo No. 18 , (May 25, 1964). 

76. Graham, E. W., "The Forces Produced by Fuel Oscillation in a Rectangular 

Tank", Douglas Aircraft Company, Inc. , Report No. SM-13748 , 

(April 1951). 

77. Lorell, J., "Forces Produced by Fuel Oscillations", Jet Propulsion 

Laboratory , Progress Report No. 20-149 , Contract No. DA-04-495-0rd 
18, (October 1951). 

78. Graham, E. W. , and Rodriguez, A. M. , "The Characteristics of Fuel 

Motion which Affect -Airplane Dynamics", Douglas Aircraft Company, 
Inc. , Report No. SM-14212 , (November 1951). 

79. Luskin, H. , and Lapin, E., "An Analytical Approach to the Fuel Sloshing 

and Buffeting Problems of Aircraft", Journal of the Aeronautical 
Sciences , (April 1952). 

80. Cooper, R. M. , "Dynamics of Liquids in Moving Containers", ARS Journal , 

(August 1960). 

81. Lawrence, H. R. , Wang, C. J., and Reddy, R. B. , "Variational Solution 

of Fuel Sloshing Modes", Jet Propulsion , (November 1958). 

82. Ring, E., Rocket Propellant and Pressurization Systems . Prentice- 

Hall, Inc., Englewood Cliffs, N. J., 1964. 


B-6 



83* Liu, F. C., "Pressure on Baffle Rings due to Fuel Sloshing in a 
Cylindrical Tank", MSFC, R-AERO-4-64, Huntsville, Alabama, 
(January 17, 1964). 

84. O'Neill, J. P., "Final Report on an Experimental Investigation of 
Sloshing", STL/TR-59 -0000 -09960 . Contract No. AF-04(647 )-309 , 
(March 4, 1960). 


B-7 



